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I.  INTRODUCTION 


The  modified  betatron1*2  and  related  devices,  like  the  stellatron3 
and  the  reversing  solenoidal  lens  (RSL)  betatron,4  are  the  subjects  of 
theoretical  and  experimental  investigations  at  present  because  of  their 
potential  as  high-current  accelerators.  For  these  devices  to  operate  suc¬ 
cessfully,  one  must  be  able  to  inject  a  high  current  beam,  in  one  or  many 
turns,  onto  a  closed  orbit  in  a  torus,  and  maintain  beam  stability  for  the 
duration  of  the  acceleration  cycle.  The  latter  may  vary  from  tens  of 
microseconds  to  several  milliseconds.  Depending  on  the  application,  it  may 
then  be  necessary  to  extract  the  beam  from  the  torus.  During  the  period  of 
our  contract  with  the  Office  Of  Naval  Research  (01  Nov  1983  -  31  Dec  1984), 
we  have  concentrated  on  the  question  of  beam  behavior  after  injection, 
i.e.,  when  a  circulating  ring  has  been  formed,  because  this  behavior  deter¬ 
mines  such  key  issues  as  the  maximum  operating  current  and  the  length  of 
the  acceleration  cycle.  One  of  the  main  concerns  in  this  regard  is  col¬ 
lective  instabilities,  including  the  negative-mass  instability,5"7  which 
can  grow  on  the  beam.  We  have  developed  a  new  analytic  theory  of  these 
instabilities  as  they  occur  in  the  modified  betatron  configuration.  Excel¬ 
lent  agreement  is  obtained  between  the  theory  and  particle  simulation 
results  over  a  wide  range  of  energy  and  wavelengths.  In  particular,  the 
short -wavelength  regime,  which  has  not  been  adequately  treated  to  date,  is 
covered  by  the  new  theory.  This  allows  us  to  obtain  the  scaling  of  growth 
rate  with  toroidal  mode-number.  Six  regimes  of  behavior  are  identified, 
according  to  the  magnitude  of  the  toroidal  magnetic  field  and  the  strength 
of  the  coupling  between  longitudinal  and  transverse  motion  (Table  1).  We 
find  that  some  subtle  geometric  effects  have  a  considerable  influence  in 
beam  stability.  Ways  to  take  advantage  of  these  effects  are  discussed.  A 
criterion  for  beam  stability  below  the  negative  mass  transition  energy  is 
given  (Eq.  65).  We  can  recover  all  the  known  limiting  cases  of  the  nega¬ 
tive  mass  dispersion  relation  from  the  theor-'.  The  model  we  use  provides 
an  analytical  framework  into  which  one  can,  at  some  future  date,  incorpor¬ 
ate  in  a  straight-forward  manner  the  effects  of  strong-focusing  fields. 
This  work  is  presented  in  Chap.  II. 
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Numerical  simulation  results  are  presented  in  Chap.  Ill  which  clearly 
illustrate  how  the  equilibrium  transverse  motion  of  the  beam  particles 
influences  beam  stability.  It  is  found  that,  for  a  given  beam  radius,  the 
most  stable  configuration  occurs  when  there  is  little  or  no  transverse 
motion  of  beam  particles.  This  is  an  important  result,  since,  as  we  shall 
show,  it  follows  that  one  does  not  necessarily  require  an  energy  spread  to 
achieve  maximum  stability.  Using  the  analytic  results  of  Chap.  II,  we  com¬ 
pute  upper  bounds  on  the  currents  stable  to  collective  instabilities. 

The  linear  theory  of  the  stellatron  and  RSL  betatron  has  not  yet  been 
worked  out.  However,  we  have  performed  several  numerical  simulations  of 
these  devices.  Results  for  the  stellatron  in  Chap.  V  show  that  in  most  of 
the  cases  looked  at,  the  quadrupole  fields  decrease,  but  do  not  eliminate, 
the  growth  of  instabilities.  We  present  some  evidence  that  the  negative- 
mass  transition  in  the  stellatron  gives  a  threshold  for  beam  stability  at 
moderate  currents  (-1  kA).  The  stability  of  the  stellatron  experiment  at 
the  University  of  California  at  Irvine8  is  addressed  in  Chap.  VI.  Limited 
results  on  the  stability  of  the  RSL  betatron  are  given  in  Chap.  VII. 

In  addition  to  looking  at  collective  instabilities,  we  have  looked  at 
some  questions  relating  to  particle  orbits  in  the  stellatron  and  RSL 
devices.  Regarding  the  stellatron,  we  find  in  Chap.  IVA  that  in  certain 
regimes  the  orbit  of  a  particle  is  sensitive  to  small  changes  in  the  ini¬ 
tial  conditions.  This  has  practical  implications  for  where  one  should 
locate  the  beam  injector  in  such  a  device.  We  have  also  done  some  self- 
consistent  (i.e.,  including  self-fields)  simulations  of  a  highly  mismatched 
beam  in  the  stellatron  (Chap.  IVB).  For  the  RSL  betatron,  we  have  looked 
at  a  system  whose  parameters  are  close  to  those  for  the  experimental  device 
under  construction  at  the  University  of  New  Mexico.  The  results  in  Chap. 
VII  show  that  by  adjusting  the  solenoidal  fields,  radial  oscillations  of 
the  beam  envelope  can  be  made  negligible  even  when  there  is  a  sizeable  gap 
between  the  solenoids. 


A  discussion  of  an  inductively-driven  instability  in  an  electron-layer 
which  is  analogous  to  that  found  in  toroidal  geometry  (Chap.  II)  is  given 
in  Appendix  A.  Appendix  B  gives  details  of  our  work  on  single-particle 
orbits  in  the  stellatron. 
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II.  NEGATIVE  MASS  INSTABILITY  LINEAR  THEORY 


In  our  last  annual  report  we  described  an  analytical  rigid-disk  model 
of  Iona  wavelength  negative  mass  instabilities  in  large  aspect  ratio  con¬ 
ventional  and  modified  betatrons  for  which  the  cavity  and  electron  beam 
minor  cross  sections  are  concentric  circles.7  This  work  represented  a  sig¬ 
nificant  advance  over  earlier  studies9.10  in  that  it  took  proper  account 
of  toroidal  corrections  to  the  electromagnetic  fields.  Agreement  between 
model  predictions  of  the  negative  mass  growth  rate  and  the  results  of  num¬ 
erical  simulations  was  quite  good.6  However,  this  agreement  showed  signs 
of  deteriorating  at  high  beam  energy  and,  of  course,  at  large  toroidal  mode 
numbers  (i.  e.,  short  wavelengths). 

During  the  past  year  we  have  relaxed  all  approximations  but  that  of  a 
rigid-disk  beam.  The  self-consistent  electromagnetic  fields  of  the  elec¬ 
tron  beam  in  the  accelerator  cavity  are  evaluated  exactly  for  arbitrary 
beam  and  cavity  dimensions  and  arbitrary  aspect  ratios.  In  addition,  the 
beam  need  not  be  centered  radially  in  the  cavity.  These  improvements  are 
made  possible  by  taking  the  beam  and  cavity  to  have  rectangular  minor  cross 
sections.  This  requirement  is  not  as  limiting  at  it  might  seem  at  first 
glance.  Three-dimensional  numerical  simulations  with  the  IVORY  code  indi¬ 
cate  that  beam  behavior  in  cavities  of  approximately  square  cross  section 
is  little  different  from  that  in  cavities  of  circular  cross  section  pro¬ 
vided  the  cross  sectional  areas  are  equal.  The  new  model  gives  even  better 
agreement  than  its  predecessor  with  negative  mass  instability  linear  growth 
rates  from  numerical  simulations.  In  particular,  discrepancies  with  the 
few  simulations  performed  at  high  energies  and  at  large  toroidal  mode 
numbers  have  been  eliminated. 

Some  new  effects  are  predicted  by  the  revised  model.  Inductive  elec¬ 
tric  fields  dominate  electrostatic  fields  at  large  toroidal  mode  numbers, 
gi/ing  rise  to  reduced  instability  growth  rates  and,  for  a  few  modes ,  even 
a  stability  window.  It  appears  that  the  stability  window  can  be  shifted  to 
low  toroidal  modes  by  moving  the  electron  beam  toward  the  inner  wall  of  the 
cavity.  We  plan  to  investigate  this  interesting  behavior  in  more  detail 
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Instabilities  in  a  modified  betatron  can  be  evaluated  in  a  manner  sim¬ 
ilar  to  that  used  for  the  conventional  betatron  provided  the  poloidal 
cyclotron  frequency,  Bq/y,  significantly  exceeds  the  toroidal  cyclotron 
frequency,  w0.  We  distinguish  between  beam  energies  comparable  to  and 
much  greater  than  the  transition  energy,  Ytr ,  defined  as  the  energy  at 
which  wz  vanishes.  Since  the  former  regime  has  no  counterpart  in  the 
conventional  betatron,  we  consider  it  later. 

For  the  modified  betatron  at  high  energies,  Eq.  (13)  reduces  to 

Q4  -  Q2  (w  2  +  z\  -  —  ?--X-a-  =  0  (52) 

V  o/ 


where 


(53) 


is  the  poloidal  drift  frequency.  (In  obtaining  these  expressions  we  have 
assumed  also  that  e  «  w02 ,  which  is  true  in  most  cases  of  interest.) 
Solutions  to  Eq.  (52)  are 


(54) 


As  before,  we  consider  the  limits  of  small  and  large  x*  For  small  x 
Eq.  (54)  reduces  to 


+  £ 


d  2/  2 
B6  h 


(55) 


When  wb2  »  e0*  Eq.  (55)  in  turn  becomes 

2  2  .2 
&  “  wb  ’  "  x''^r 


(56) 
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For  small  x»  i.  e.,  weak  coupling  between  longitudinal  and  transverse 
beam  modes, 

*  ( T -T?)2  7  (%  -  H  ve2)  «  1  <48> 

r 

the  two  pairs  of  roots  in  Eq.  (47)  reduce  to 

2  2 
£  *  “>r  .  -  x/^r 

The  former  represents  transverse  (in  r)  oscillations,  which  are  stable,  and 
the  latter  represents  predominantly  longitudinal  oscillations,  one  of  which 
is  unstable  with  a  growth  rate 

r  ■  I  [rhr  7  (93  •  94  v(.2)]1/2  <5»> 

if  x  (equivalently,  e0)  is  positive.  Eq.  (50)  becomes  the  usual  negative 
mass  instability  growth  rate  expression  when  g3  =  g^.14,15  However,  the 
instability  vanishes,  if  x  is  negative.  This  stabi 1 ization  of  the  negative 

mass  effect  is  caused  by  inductive  coupling  of  the  beam  to  the  accelerator 

cavity,16.17  which  is  discussed  further  in  the  next  section. 

With  increasing  i  the  inequality  (48)  eventually  is  reversed,  giving 
rise  to  strong  coupling  between  the  longitudinal  and  transverse  modes.  One 
unstable  hybrid  mode  with  a  growth  rate 

a172  rv  /  2\i174 

r  ■  V  [7  (93  -  94  ve  )]  <51> 

occurs  if  x  is  positive,  and  two  unstable  modes  with  growth  rates  reduced 

from  Eq.  (51)  by  /?  occur  otherwise.  Numerical  calculations  performed  to 
date  suggest  that  the  latter  possibility  is  much  more  common.  Note  that 
the  growth  rate  at  large  i,  Eq.  (51),  scales  as  A1/2,  in  contrast  to  the 
growth  rate  at  small  SL ,  Eq.  (50),  which  scales  as  i. 


(49) 


< 
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Although  the  precise  values  of  g^  and  Q2  are  not  critical,  the 
three  remaining  coefficients  must  be  determined  quite  accurately.  Expand¬ 
ing  Eq.  (40)  and  dropping  small  terms  proportional  to  Q2  and  v2  yields 

x  "  (g5  "  g4  ve2)  +  ^  7  (g3  ‘  g4  ve2)  (43) 


Thus,  the  coupling  term  x  between  longitudinal  and  transverse  beam  modes 
depends  sensitively  upon  the  differences  between  g5  and  g^,  and  g3  and  g^. 
Numerical  evaluations  of  these  factors  indicate  that  the  second  term  in 
Eq.  (43)  is  more  often  the  larger,  in  which  case 


x 


7?  eo 


(44) 


where  e0  is  e  evaluated  at  Q  =  0. 

Let  us  now  specialize  to  a  high  current  conventional  betatron,  for 
which  the  general  dispersion  relation  is 

(q2  -  wr2)  (g2  -  e)  =  x  (45) 

For  an  equilibrium  to  exist  in  the  conventional  betatron  at  high  current,  y 
must  be  sufficiently  large  that  the  first  term  in  Eq.  (38)  dominates  the 
second.  Then,  Eq.  (45)  becomes  approximately 

G4  -  oor2C2  -  x  =  0  (46) 

with  solutions 
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less  accurate  but  still  tolerable.  This  is  not  true  of  Eq.  (17),  however, 
due  to  a  strong  cancellation  between  electrostatic  and  inductive  contribu¬ 
tions.  With  these  approximations  the  terms  occurring  in  the  dispersion 
relation  take  the  form 
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(40) 
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where  u  is  the  beam  current  normalized  to  17  kA.  For  a  square  beam  roughly 
centered  in  a  square  cavity,  the  five  geometrical  factors  are  approximately 

gx  -  g2  ■  1  (41) 

93  "  g4  "  %  *  7  +*n(ir)  (42) 

as  can  be  verified  numerically.  Ac  is  the  cavity  cross  sectional  area, 
defined  similarly  to  A. 


Recall  that  a,  which  typically  is  of  order  unity,  occurs  only  in  Eq.  (15). 
The  eigenfunctions  entering  into  Eq.  ( 29 ) - ( 32)  are 


X 

9E  ‘ 

ce1/2  [J*(V)  ¥ji(v?)  -  J*(V2)  Y.<(V)] 

(33) 

X 

%  = 

ci1/2  [j*(V)  v;(V2)  -  j;(V2)  v*(V)] 

(34) 

The  quantities  J*  and  Y*  are  standard  Bessel  functions  of  the  first  and 
second  kinds.12  A  prime  indicates  differentiation  with  respect  to  the 
argument.  The  normalization  constants  in  these  expressions  are13 

CE  ■  1  [r2  teVE  <35) 

CB  ■  H(r  - 

Accurate  solutions  of  Eqs.  (13)-(36)  clearly  must  be  obtained  numeric¬ 
ally.  However,  a  qualitative  picture  of  instability  scaling  can  be  devel¬ 
oped  with  surprising  ease. 

B.  Approximate  Scalings  of  Instability  Growth  Rates 

The  negative  mass  and  related  instabilities  typically  occur  at  fre¬ 
quencies  approximated  by  Q  »  0.  Thus,  provided  is  well  below  the  elec¬ 
tromagnetic  cutoff  for  toroidal  mode  A  of  the  betatron  cavity,  which  is 
usually  the  case,  we  can  safely  replace  u>  by  in  the  denominators  of  the 
field  contributions  (the  sums)  in  Eqs.  (14)-(27).  Making  the  same  substi¬ 
tution  in  the  numerators  of  the  field  contributions  in  Eqs.  ( 14 )- ( 16)  is 
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r2  -  ri 


n  =  0 


n  >  0 


(The  lowest  electric  eigenvalue  is,  of  course,  not  included  in  the  sums.) 
However,  as  i  increases  and  the  magnitude  of  the  lowest  magnetic  eigenvalue 
approaches  that  of  the  second  lowest,  the  values  of  the  low-lying  electric 
and  magnetic  eigenvalues  are  pushed  up. 

The  radial  form  factors  are  expressed  as  integrals  over  the  vacuum 
eigenfunctions,  because  simpler  expressions  do  not,  in  general,  exist. 
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with,  from  Eq.  (6), 


P  =  Y  <*> 


o  -  VJt 


E  Kl 


k  \2  \  *  \ 

PE  PE 


n  ,  o.  2 

k,°XE  k  XE 


(19) 


Sums  in  the  preceding  six  equations  are  performed  over  the  axial  and 
radial  eigenvalues  of  the  vacuum  electromagnetic  fields  in  the  toroidal 
cavity.  The  axial  eigenvalues  are  given  by 


.  2u 
k  =  —  n 

a 


n  =  [0,  ®) 


(20) 


k  k 

The  axial  form  factors  pE  and  Pg  are  obtained  by  integrating  the  corres¬ 
ponding  eigenfunctions,  sines  and  cosines,  over  the  beam  cross  section. 
r 
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0 

n  even 

(22) 

eigenvalues  satisfy  the  usual 

vacuum  dispersion 

relations: 

:  rl)  Y*(  XE  r2)  " 

Ja(  ^ Er 2) 

v*(Vi)  -  0 

(23) 

Brl)  V(  XBrz)  ' 

'  (A^Br2) 

V(Vl)  ’  0 

(24) 

For  low  mode  numbers,  l ,  these  eigenvalues  can  be  approximated  as 

0 

2n 


f2  '  rl 


n 


n  =  0 

n  >  0 


(25) 


1R 


the  longitudinal  oscillation  frequency  (squared)  by 


and  the  coupling  coefficient  by 


Y 


B°  is  the  applied  vertical  magnetic  field,  and  n  is  the  betatron  index 
The  vertical  field  is  determined  from  radial  force  balance,  Eq.  (1). 


For  the  calculations  outlined  below,  it  is  convenient  to  express 
Maxwell's  equations  for  the  perturbed  field  quantities  as  a  vector  wave 
equation  for  the  electric  fields. 


6E  =  -  6J  -  7  6p 
ot 


Once  Eq.  (11)  is  solved,  the  magnetic  fields  are  determined  simply  from 


It  6B 


-  Vx6E 


The  solution  of  Eq.  ( 1 )-( 12)  is  straightforward  but  very  time-consum¬ 
ing.  It  proceeds  as  follows.  The  perturbed  quantities  are  taken  to  have 
azimuthal  and  temporal  variations  of  the  usual  form  exp[i(£0  -  wt)].  Eq. 
(11)  together  with  its  source  terms,  Eq.  (7)-(10),  are  then  solved  for  the 
perturbed  electric  fields  using  vector  Green's  functions  in  cylindrical  co¬ 
ordinates.11  Next,  the  perturbed  magnetic  fields  are  obtained  by  direct 
substitution  in  Eq.  (12).  The  equilibrium  fields  are  determined  in  a  simi¬ 
lar  manner.  (Note  that  the  electric  and  magnetic  fields  also  can  be  der¬ 
ived  starting  from  scalar  Green's  functions  for  Ez  and  Bz.  However,  we 
have  found  the  resulting  infinite  series  to  be  more  prone  to  numerical  dif¬ 
ficulties.)  Eliminating  the  various  field  components  from  Eq.  ( 3 ) - ( 5 )  then 
leaves  a  Hermitian  matrix  equation  for  the  perturbed  beam  centroid  loca¬ 
tion.  Its  determinant  is  the  desired  dispersion  relation. 

The  dispersion  relation,  cast  in  the  form  used  in  our  earlier 
work6  »7 ,  is 


(°2  -  “/)  {f  -  “r2  - 


^)-2  V  ■  » 


where  Q  =  oj  -  lu  and  w  =  VL/R.  The  axial  and  radial  oscillation  frequen- 
O  0  w 

cies  are  given  by 
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(6) 


Perturbed  quantities  are  preceded  by  a  delta  (e.g.,  6z),  while  unperturbed 
quantities  are  not.  Total  time  derivatives  of  perturbed  beam  quantities 
are  represented  by  dots  above  those  quantities  (e.g.,  6z).  B0  is  the  tor¬ 
oidal  magnetic  field,  averaged  over  the  beam  cross  section. 


Perturbed  beam  current  and  charge  densities  arising  from  these  motions 
are 


6J2  - 

p  6z 

(7) 

6Jp  = 

• 

p  6r 

(8) 

7  6Je 

• 

V 

=  p  66 

+  IT  6p 

(9) 

and 

6p  *  '  I?  p  6z  '  F  t?  r  p6r  "fe  p6e  <10> 

With  the  exception  of  Eq.  (9),  which  has  been  modified  to  take  account  of 
the  beam  azimuthal  velocity  variation  with  radius,  these  equation  are  taken 
directly  from  Ref.  6.  It  is  easy  to  show  that  they  satisfy  the  continuity 
equation. 
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balance  equation  is  given  by 


J  /  (Er  +  r  r  Bz)  dA  +  *  V  -  0  <*> 

In  arriving  at  Eq.  (1),  we  have  taken  the  beam  density  to  be  constant  and, 
consistent  with  the  rigid-disk  approximation,  the  beam  velocity  to  vary 
linearly  with  radius.  R  is  the  charge-weighted  equilibrium  radius,  Ve  is 
the  velocity  at  that  location,  A  is  the  cross  sectional  area  of  the  beam, 
and  dA  =  r  dr  dz. 

R  -  7  f 

Note  that  the  speed  of  light  and  the  electron  charge  and  mass  are  scaled 
out  of  our  equations.  This  has  the  effect  of  measuring  distance  and  time 
relative  to  an  arbitrary  scale  length,  and  potentials  relative  to  the  elec¬ 
tron  rest  energy. 

The  linearized  equations  for  beam  motion  about  its  equilibrium  posi¬ 
tion  are  obtained  in  a  similar  manner.7.9 
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(including  IVORY  simulations)  in  the  coming  year.  (The  analogous,  but  much 
simpler,  case  of  inductive  effects  in  a  fast  rotating  e-layer  immersed  in  a 
strong  azimuthal  magnetic  field  is  treated  briefly  in  Appendix  A.  Besides 
providing  a  check  on  the  betatron  results,  it  permits  more  insight  into  the 
underlying  physics.)  High-frequency  electromagnetic  effects  on  betatron 
stability  also  were  investigated  but  were  found  to  be  small  for  cavity 
cross  sections  of  moderate  size. 

The  remainder  of  this  chapter  is  organized  as  follows.  In  Section  A 
the  beam  equations  of  motion  are  specified,  and  the  electromagnetic  fields 
are  evaluated  by  a  Green's  function  expansion.  Although  the  resulting 
dispersion  relation  is  algebraically  very  complicated,  useful  analytical 
growth  rate  expressions  can  be  derived  in  several  limiting  cases,  described 
in  Section  B.  The  numerical  results  in  Section  C  support  and  extend  our 
analytic  findings  for  1  and  10  kA  beams  over  a  wide  range  of  magnetic 
fields  and  mode  numbers.  IVORY  comparisons  also  are  provided.  Our  results 
are  summarized,  their  implications  discussed,  and  directions  of  further 
linear  stability  research  suggested  in  Section  D. 

A.  Derivation  of  the  Linear  Dispersion  Relation 

We  consider  a  rectangular  beam  in  a  rectangular  torus  with  dimensions 
as  shown  in  Fig.  1.  In  general,  the  beam  cross  section  need  not  be  the 
same  shape  as  that  of  the  torus,  and  the  beam  need  not  be  centered  radi¬ 
ally.  Axial  centering  is,  of  course,  required  for  a  static  equilibrium. 

As  in  earlier  works,6  »7»9*  10  we  treat  the  beam  as  a  string  of  non- 
rotating  rigid  disks.  Thus,  the  internal  dynamics  of  the  beam,  and  temper¬ 
ature  effects  in  particular,  are  ignored.  (Negative  mass  instability  sta¬ 
bilization  by  an  electron  energy  spread  is  discussed  in  Chapter  III  of  this 
report,  and  el sewhere.b »7 )  Within  this  context  the  beam  centroid  equa¬ 
tions  of  motion  are  obtained  by  integrating  the  single  particle  equations 
over  the  beam  cross  section.  Thus,  the  beam  equilibrium  radial  force 
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Note  that  the  second  root  in  Eq.  (56)  is  identical  to  that  in  Eq.  (49). 
Consequently,  the  toroidal  magnetic  field  has  no  effect  on  the  negative 
mass  instability  whenever 


l2 
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r5 


2  2, 
10  / 
:  r 


(s3  -  =4  ve2) 


(57) 


Very  large  currents  and  applied  fields  are  needed  to  violate  this  inequal¬ 
ity,  except  at  small  y. 


In  the  opposite  limit  Eq.  (55)  simplifies  to 
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*  e 


U2  X/£0 
°*  Be2^2 


(58) 


The  second  of  these  roots  can  be  simplified  by  means  of  Eq.  (44)  to  the 
usual  negative  mass  instability  growth  rate  expression  for  very  large 
toroidal  fields,18 

r  «  (1  -  n)"1^  y  <*>g  (59) 


which  is  independent  of  both  A  and  v.  The  first  root  in  Eq.  (58)  is  purely 
longitudinal.  Instability  results  for  e0  <  0,  i.  e.,  for  inductive  elec¬ 
tric  fields  dominating  over  electrostatic. 

Returning  to  Eq.  (54),  we  find  that  the  large  x  limit  corresponds  to  a 
strong  coupling  between  longitudinal  and  transverse  modes,  just  as  in  the 
conventional  betatron.  The  growth  rate  here  is  simply  that  described  by 
Eq.  (51),  reduced  by  of  order  ( Bz/B0  J1 /2 . 
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The  various  cases  we  have  considered  thus  far  for  both  the  conven¬ 
tional  and  the  modified  betatron  are  categorized  in  Table  1  according  to 
the  magnitudes  of  u>b  and  x. 


It  is  more  difficult  to  obtain  useful  instability  growth  rate  esti¬ 
mates  for  the  modified  betatron  at  low  beam  energies,  because  the  ordering 
of  most  terms  in  the  dispersion  relation  is  not  uniform.  Replacing  e  by 
co  often  is  particularly  dubious.  Nonetheless,  some  progress  can  be  made 
for  u)22  near  zero,  provided  wr2  also  is  small  there.  In  this  case  the 
dispersion  relation  reduces  approximately  to 


Q4  -  Q2  e 
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which  has  as  a  formal  solution 
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As  u>z  goes  to  zero,  so  does  one  root  of  Eq.  (62),  and  replacing  e 


by  e0  becomes  permissible. 
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TABLE  1.  Approximate  solutions  to  the  dispersion  relation,  Eq.  (13),  for 
conventional  betatrons  and  for  modified  betatrons  well  above  the 
transition  energy.  Solutions  are  categorized  by  the  toroidal 
field  strength  B0  and  the  coupling  coefficient  x • 
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Eq.  (63)  is  similar  to  the  strong  magnetic  field  expression,  Eq.  (58),  but 
with  a  correction  term  in  the  denominator.  When  the  denominator  is 
positive,  Fq.  (63)  predicts  instability  above  the  transition  energy  (uz2 
>  0)  and  stability  below.  Note,  however,  that  further  below  the  transition 
energy  the  second  term  in  the  square  root  in  Eq.  (62)  may  dominate  the 
first,  i.  e.. 


causing  the  two  unstable  modes  often  found  numerically  at  low  energies 
(Refs.  6,  7,  9,  10).  A  large  enough  magnetic  field  eliminates  the 
instability  there. 

Be  »  Ytr2  w0  (65) 

Incidentally,  such  a  small  toroidal  magnetic  field  that  the  denomin¬ 
ator  of  Eq.  (63)  becomes  negative  causes  the  root  described  by  this  equa¬ 
tion  to  be  stable  above  the  transition  energy  but  unstable  below.  The 
other  root  in  Eq.  (62)  is  unstable  throughout. 

C.  Numerical  Solutions  of  the  Linear  Dispersion  Relation 

A  computer  program,  BTRSQ,  has  been  written  to  evaluate  numerically 
the  roots  of  the  dispersion  relation.  Sums  appearing  in  Eq.  (13)  are  per¬ 
formed  over  the  first  twenty  to  forty  axial  and  radial  eigenvalues  and  gen¬ 
erally  converge  rapidly.  Good  accuracy  for  u>z2  and  wr2,  however,  is 
achieved  for  a  moderate  number  of  terms  only  when  the  numbers  of  axial  and 
radial  eigenvalues  used  are  approximate  integer  multiples  of  twice  the 
cavity-to-beam  ratios  of  the  axial  and  radial  cross  sectional  dimensions. 


respectively.  The  roots  themselves  of  the  dispersion  relation  are  deter¬ 
mined  by  Muller's  method.19  Typically,  five  seconds  of  CDC-7600  computer 
time  is  needed  to  initialize  an  equilibrium  configuration,  after  which 
roots  can  be  obtained  at  the  rate  of  about  one  per  second. 

Calculations  were,  for  the  most  part,  performed  using  the  equilibrium 
parameters  listed  in  Table  2  and  toroidal  mode  numbers  in  the  range  one  to 
forty.  The  parameters  were  chosen  to  bracket  those  of  the  modified  beta¬ 
tron  under  development  at  the  Naval  Research  Laboratory.20 

Preliminary  insight  can  be  developed  by  observing  the  scaling  with 
beam  and  cavity  dimensions  and  with  toroidal  mode  number  l  of  the  geomet¬ 
rical  factors  g^  through  g,.  introduced  by  Eqs.  (37)  -  (40).  For 

simplicity,  we  eliminate  u>  from  the  sums  in  Eqs.  (14)  -  (17),  replacing  it 
by  Jt/R  in  the  denominators  and  by  Xco0  in  the  numerators  of  the  terms  in 
the  sums.  The  resulting  expressions  for  g^,  g4,  and  g^  are 

approximately  independent  of  beam  energy,  and  the  energy  dependence  of 
g^  and  g£  can  be  made  explicit  as 

«i  ■  1,2  (V-  ^V)  (66> 

g2  =  y2  (g2a-  92\2)  (67) 

Note  that  in  deriving  Eq.  (67)  we  have  dropped  an  extremely  small  term  pro¬ 
portional  to  Er2 ,  because  it  is  inconsistent  with  the  scaling  shown. 

Table  3  lists  geometrical  factors  for  a  range  of  Jl  values,  based  on 
the  configuration  in  Table  2.  Their  magnitudes  are  approximately  as 
expected,  namely  unity  in  the  first  four  columns,  and  Eq.  (42),  here  3.72, 
in  the  last  three  columns.  More  importantly,  relative  differences  among 
related  factors  are  small.  The  few  percent  difference  between  g^a  and  g^, 
for  instance,  indicates  that  is  essentially  independent  of  y  for  y  < 
10.  At  much  higher  energies  g^  increases  as  y2 ,  but  at  such  energies 


TABLE  2.  Typical  betatron  parameters  used 

in 

numerical  solutions  of  the 

dispersion  relation,  Eq.  (13). 

Beam  Current 

1,  10  kA 

Beam  Energy 

1-30  MeV 

Toroidal  Magnetic  Field 

0, 

1,  10,  100  kG 

Index  (n) 

0.5 

Beam  Inner  Radius  (w^) 

98.24  cm 

Beam  Outer  Radius  (W2) 

101.76  cm 

Beam  Axial  Half-Width  (b) 

1.76  cm 

Torus  Inner  Radius  (ri) 

91.2  cm 

Torus  Outer  Radius  (r2) 

108.8  cm 

Torus  Axial  Half-Width  (a) 

8.8  cm 

TABLE  3. 

Approximate  geometrical  factors  computed  for 
Table  2  at  various  toroidal  mode  numbers. 

the  parameters  in 

Jl 

„  a 
gl 

glb 

„  a 
g2 

g2b 

g3 

g4 

g5 

1 

1.094 

1.090 

1.094 

1.093 

3.543 

3.538 

3.571 

5 

1.094 

1.090 

1.095 

1.093 

3.542 

3.539 

3.570 

10 

1.095 

1.090 

1.097 

1.091 

3.542 

3.541 

3.568 

15 

1.095 

1.089 

1.100 

1.089 

3.541 

3.545 

3.567 

20 

1.096 

1.088 

1.103 

1.084 

3.540 

3.551 

3.572 

25 

1.097 

1.086 

1.104 

1.075 

3.540 

3.559 

3.589 

30 

1.098 

1.084 

1.104 

1.062 

3.540 

3.570 

3.627 

35 

1.098 

1.080 

1.098 

1.041 

3.542 

3.586 

3.700 

40 

1.097 

1.075 

1.086 

1.009 

3.546 

3.607 

3.828 

30 


self-fields  have  no  effect  on  wz2.  Likewise,  the  small  differences 
between  and  g^  justifies  our  earlier  approximation  that  the  term 
in  Eq.  (43)  linear  in  Q  usually  can  be  ignored. 

Differences  between  g^  and  g^,  on  the  other  hand,  have  a 
strong  impact  on  the  beam  stability,  because  they  can  affect  the  magnitude 
and  even  the  sign  of  e  and  x*  At  small  A,  g^  slightly  exceed*  g^, 
and  as  a  consequence  negative  mass  growth  rates  fall  off  at  high  energies 
as  y'1/2  rather  than  the  usually  cited  y~3/2.  The  difference  between 
the  two  geometry  factors  decreases  with  increasing  toroidal  mode  number,  so 
that  at  A  =  11  they  are  equal.  At  higher  A,  g^  exceeds  g^,  stabil¬ 
izing  the  beam  at  sufficiently  large  energies.  As  we  shall  soon  see,  this 
reprieve  from  the  negative  mass  instability  is  short-lived.  For  still 
larger  A  the  beam  encounters  the  large  x  regime  indicated  by  the  second, 
fourth,  and  sixth  cases  in  Table  1,  and  the  hybrid  instability  discussed  in 
the  preceding  section  occurs. 

The  coefficients  g^  and  g^  represent  the  electrostatic  and 
inductive  electric  self-field  contributions  to  the  beam  longitudinal 
dielectric  function.  Electrostatic  forces  always  exceed  inductive  in  a 
straight  drift  tube,  if  no  slow-wave  structure  is  present.  Evidently,  the 
larger  inductive  field  in  a  betatron  is  due  to  curvature  of  the  beam  and 
cavity.  Ordinarily,  one  would  expect  a  dominant  inductive  field  to  cause 
unstable  longitudinal  bunching,  as  in  a  Cherenkov  maser.21  The  reverse, 
however,  is  true  of  a  beam  particle  with  effectively  negative  mass,  for 
which  excess  electrostatic  fields  cause  instability.  To  corroborate  this 
picture,  we  have  compared  the  stability  of  beams  in  betatrons  to  that  of 
rotating  electron  layers,  a  simpler  problem.  Both  were  immersed  in  very 
strong  azimuthal  guide  fields  to  suppress  negative  mass  effects.  Inductive 
fields  were  found  to  dominate  electrostatic  in  both  cases  at  high  A 
values.  Behavior  of  the  electron  layer  is  described  more  fully  in  Appendix 


Table  4  illustrates  the  changes  in  the  geometrical  coefficients  for  a 
=  1  as  the  position  of  the  beam  centroid  is  varied.  When  the  beam  is  near 
either  the  inner  or  tne  outer  wall  of  the  cavity,  g1  is  much  reduced 
due  to  shorting  of  the  axial  self-field  restoring  forces,  and  g?  is 
correspondingly  increased. 

TABLE  4.  Approximate  A  *  1  geometrical  factors  computed  for  the  parameters 
in  Table  2  but  with  the  beam  offset  radially  by  various  amounts. 


R(cm) 

9l3 

b 

*1 

s2a 

b 

92 

93 

94 

95 

107 

0.136 

0.134 

12.768 

13.282 

1.567 

1.557 

1.326 

103 

0.851 

0.847 

1.796 

1.875 

3.326 

3.319 

3.357 

100 

1.094 

1.090 

1.094 

1.093 

3.543 

3.538 

3.571 

97 

0.905 

0.903 

1.901 

1.810 

3.326 

3.324 

3.314 

95.67 

0.684 

0.684 

3.133 

2.962 

3.037 

3.037 

2.986 

93 

0.158 

0.159 

14.601 

13.936 

1.577 

1.587 

1.260 

Hence,  beam  stability  at  low  energies  can  be  affected  strongly  by  the  cen¬ 
tering  of  the  beam.  In  addition,  g^  grows  relative  to  g3  as  R  is 
decreased.  The  two  become  equal  at  R  =  95.67  cm,  and  for  smaller  radii  the 
inductive  field  dominates.  This  suggests  that  improved  stability  at  low  A 
can  be  achieved  for  a  high  energy  beam  by  placing  it  near  the  inner  wall  of 
the  torus. 

Figures  2-5  depict  A  =  1,  10,  20,  and  30  negative  mass  instability 
growth  rates  for  10  kA  and  1  kA  beams  in  1  kG  guide  fields,  determined 
directly  from  Eq.  (13).  Results  for  A  =  1  and  10  with  y  >  10  agree  well 
with  Eq.  (50).  The  growth  rate  scales  linearly  with  A  and  /v.  As  predicted 
above,  the  growth  rates  at  high  energies  vary  roughly  as  y1/2  for  A  =  1 
and  as  y3/2  for  A  =  10.  Scaling  of  growth  rates  at  low  energies  is  not 
so  clean  but  appears  to  go  as  (A/v)1/2.  Although  the  two  unstable  modes 
below  ytr  appear  to  join  directly  onto  the  one  unstable  mode  above  ytr» 
a  closer  examination  indicates  that  a  tiny  gap  exists,  consistent  with  the 
discussion  accompanying  Eq.  (63). 
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Figure  3.  Growth  rates  of  the  £  -  10  negative  mass  instability  versus 
energy  for  10  kA  and  1  kA  beams  in  1  kG  guide  fields. 
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When  l  >  11  the  standard  negative  mass  instability  is  cut  off  by 
inductive  effects  for  %  negative,  i.  e.,  y  greater  than 

''co  •  (»  -  83/94)'1/2  <68> 

The  cutoff  shown  for  the  1  kA  beam  at  X  =  20  in  Fig.  4  satisfies  Eq.  (68) 
well.  Note  also  the  first  appearance  of  the  hybrid  negative  mass  instabil¬ 
ity  in  a  narrow  band  around  y  =  24.  The  hybrid  mode  is  much  stronger  in 
the  10  kA  beam  and  merges  directly  into  the  usual  branch  of  the  instability 
at  Yco •  The  growth  rates  are  described  accurately  by  Eq.  (54),  but  not 
yet  by  Eq.  (60)  because  x  is  not  sufficiently  large.  In  Fig.  5  for  Z  =  30, 
yc0  «  12;  the  hybrid  negative  mass  instability  is  well  developed  for  both 
10  kA  and  1  kA. 

Figure  6  is  intended  to  illustrate  the  mode  couplings  that  give  rise 
to  the  different  branches  of  the  negative  mass  instability.  It  presents 
the  Doppler-shifted  real  frequencies  Q  corresponding  to  the  growth  rates  of 
the  10  kA  beam  in  Fig.  4.  Below  Ytr  and  above  Yco  the  longitudinal  and 
transverse  modes  are  strongly  coupled,  and  two  unstable  branches  with  com¬ 
parable  growth  rates  typically  exist.  The  standard  negative  mass  instabil¬ 
ity  occurs  within  the  intervening  energy  range,  resulting  from  a  non¬ 
resonant  distortion  of  the  longitudinal  modes  by  curvature  effects. 

Instability  growth  rate  scaling  with  Be  is  depicted  in  Fig.  7 
and  8  for  the  10  kA  beam  at  Z  =  1  and  20.  Results  for  1  kG,  10  kG,  and  100 
kG  guide  fields  are  shown.  Corresponding  data  for  the  conventional  beta¬ 
tron  are  not  given  for  Z  =  1,  and  for  Z  -  20  below  yco.  because  they  are 
indistinguishable  from  the  1  kG  results  at  energies  for  which  an  equilib¬ 
rium  exists  (above  about  13  MeV).  The  conventional  betatron  is  stable 
above  Yco  for  *  =  20,  because  inequality  (48)  is  not  yet  strongly 
reversed.  In  general,  we  expect  the  guide  field  to  have  negligible  influ¬ 
ence  on  beam  stability  for  low  A  (e.  g..  Fig.  7)  so  long  as  Eq.  (57)  is 
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Figure  14.  Cold  beam  gri 
1  kA  beam  ve 


It  is  important  to  note  that  the  upper  bound  on  the  frequency  spread 
given  by  case  (a)  applies  regardless  of  the  type  of  equilibrium  with  which 
one  is  dealing.  There  are  various  ways  of  trying  to  achieve  this  upper 
bound.  One  way,  as  we  have  just  pointed  out,  is  to  have  a  Finn-Manheimer 
equilibrium  with  sufficiently  small  wpr.  Alternatively,  one  can  use  the 
RSL  configuration,  which  as  we  shall  see  in  Chap.  VII,  gives  a  beam  which 
has  maximum  frequency  spread  at  all  energies.  Finally,  one  may  introduce 
an  energy  spread  on  the  beam.5.9  However,  self-consistent  modified  beta¬ 
tron  equilibria  which  include  an  energy  spread  have  not  been  adequately 
studied  to  date.  Limited  simulations7  indicate  that  it  may  be  difficult  to 
achieve  an  equilibrium  whose  radius  is  determined  by  energy  spread  rather 
than  by  betatron  oscillations.  Analytic  calculations  of  beam  stability  in 
the  presence  of  an  energy  spread  generally  assume  that  the  beam  radius  is 
determined  solely  by  the  energy  spread.5  This  is  a  very  optimistic  assump¬ 
tion. 

We  now  look  at  some  IVORY  simulations  which  illustrate  the  effect  of 
frequency  spread  on  beam  stability.  For  all  the  examples  considered,  A0 
■»  8*10-5  cm-1,  since  a  =  2  cm,  R  =  100  cm.  First,  we  look  at  a  1  kA  beam 
at  y  =  6  in  a  5  kG  toroidal  field,  other  parameters  being  given  in  Table 
5.  The  diamagnetic-paramagnetic  transition  for  this  beam  occurs  approxi¬ 
mately  at  y  =  8.  We  simulated  the  i  =  20  mode  on  this  beam  using  IVORY. 
The  cold-beam  growth  rate  for  this  mode,  from  Fig.  14,  is  6*10-4  cm-1, 
while  the  simulation  growth  was  much  smaller,  at  <_  1*10"5  cm-1.  From  a 
computer  movie  of  the  simulation,  the  poloidal  rotation  frequency  was  meas¬ 
ured  to  be  5X10"4  cm-1,  which  is  significantly  smaller  than  JlA0  c 
1.6X10-3  cm-1.  One  would  therefore  expect  case  a)  above  to  apply,  so  that 
the  growth  rate  is  reduced  by  the  amount  £A0,  which  is  more  than  enough 
to  stabilize  the  mode,  as  observed.  The  Jt  =  1  mode  on  this  beam  is  pre¬ 
dicted  to  be  unstable,  however.  The  cold  beam  growth  rate  from  Fig.  14  is 
1.8X10~4  cm-1,  which  is  greater  than  A0.  We  have  not  verified  with  IVORY 
that  the  beam  is,  in  fact,  unstable,  because  the  code  would  have  to  be  run 
a  rather  long  time  to  give  a  good  estimate  of  the  small  expected  growth 
rate,  1.0x10-4  cm-1.  As  an  aside  which  is  relevant  to  experimental  efforts 
going  on  at  NRL  and  elsewhere,  we  note  from  Fig.  14  that  the  growth  rates 


rate  of  the  instability  on  this  beam  is  accurately  given  by  the  cold-beam 
growth  rate,  as  we  shall  see  below.)  In  Fig.  13(c),  y  =  20,  so  that  the 
beam  is  near  the  transition,  and  therefore  is  rotating  more  slowly.  In 
this  case,  we  see  that  the  spreading-out  of  the  bunches  in  the  toroidal 
direction  is  much  larger.  This  figure  shows  clearly  that  the  spreading-out 
occurs  because  particles  at  smaller  radius,  with  a  shorter  distance  to 
travel,  move  ahead  relative  to  those  at  larger  radius. 

Let  us  now  proceed  with  the  calculation  of  an  effective  value  of  a  for 
the  two  cases  mentioned.  For  case  a),  a  particle  at  radius  r  =  R  +  6r  has 
circulation  frequency  Q=u0( l+6r/R) .  Assuming  that  the  beam  has  uniform 
density  profile,  the  distribution  of  frequencies  is  then  fb  *  ((Q-w0)2 
-  w02a2/R2  )l  /2 ,  where  a  is  the  beam  radius.  Unlike  f|_  above,  this 
distribution  has  a  sharp  cutoff,  resulting  from  the  finite  extent  of  the 
beam.  Nevertheless,  we  can  compute  an  effective  A  for  the  beam  by,  for 
example,  demanding  that  the  half-widths  of  the  cumulative  distributions  be 
the  same  for  fj_  and  fb-  In  physical  space,  this  means  that  one  is 
specifying  a  radial  interval  R  -  br^  to  R  +  Ar^  and  requiring 
that  half  the  beam  particles  be  within  this  band.  This  criterion  leads  to 
the  relation  A  *>  0.4ojQa/R  =  A0.  If  one  applies  a  similar  criterion  to 
compute  an  effective  width  for  case  b),  one  finds  that  the  cumulative  dis¬ 
tribution  corresponding  to  the  Lorentzian  is  divergent,  i.e.,  a  finite  beam 
width  cannot  be  defined  for  fi_.  We  will  not  dwell  on  this  difficulty 
here  but  simply  note  that  for  those  simulations  we  have  performed  where 

case  b)  applies,  the  frequency  spread  present  has  a  negligible  effect  on 

the  growth  rate.  Presumably  the  effect  of  the  small  frequency  spread  would 
become  important  only  at  X  values  significantly  higher  than  those  we  have 
considered  to  date  (viz.  X  =  1-20). 

For  the  Finn-Manheimer  equilibria,  the  beam  has  some  poloidal  rotation 
unless  one  is  right  at  the  transition  energy.  It  is  therefore  reasonable 

to  ask  how  we  decide  whether  a  particular  simulation  run  falls  under  case 

(a)  or  case  (b)  above.  It  is  easy  to  show  that  if  JIA0  >  wpr  where 
wpp  is  the  poloidal  rotation  frequency,  and  X  is  the  mode  number  we  are 
simulating,  then  the  beam  falls  under  case  (a).  In  the  opposite  limit, 
case  (b)  applies. 
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Illustration  of  dependence  of  toroidal  rotation  frequency- 
spread  on  beam  energy.  Initial  particle  positions  are  shown 
in  (a).  The  beam  in  (c)  is  close  to  the  transition  energy 
=  18,  so  the  toroidal  mixing  is  faster  than  in  (b). 


III.  UPPER  BOUNDS  ON  STABLE  CURRENTS  IN  THE  MODIFIED  BETATRON 

In  the  analytic  model  of  Chap.  II,  all  particles  have  the  same  circu¬ 
lation  frequency  about  the  major  axis  since  the  beam  is  assumed  to  be  made 
up  of  rigid  disks.  One  would  expect  a  spread  in  circulation  frequencies  of 
the  particles  to  lead  to  a  reduction  in  growth  rates  through  Landau  damp- 
ing.  *  If  one  assumes  that  the  frequency  distribution  is  Lorentzian, 
f L (Q )  a  l/((c-wo)2  +  a2)»  where  w0  is  the  average  (matched)  circula¬ 
tion  frequency,  and  a  is  a  constant,  then  a  simple  calculation5  shows  that 
the  growth  rate  of  an  instability  with  toroidal  mode-number  H  should  be 
reduced  by  JlA.  A  Lorentzian  is  a  somewhat  unrealistic  distribution,  having 
long  tails,  but  more  realistic  distributions  are  much  more  difficult  to 
treat  analytically.  The  Lorentzian  is,  therefore,  used  to  give  rough  esti¬ 
mates.  We  now  look  at  how  one  can  make  a  reasonable  choice  of  the  value  of 
A  for  the  types  of  equilibria  we  have  been  simulating  with  IVORY.  Since  we 
are  concerned  for  the  most  part  with  relativistic  beams,  all  the  particles 
have  about  the  same  speed.  A  spread  in  circulation  frequencies  can  there¬ 
fore  result  only  from  a  spread  in  path  lengths  of  the  particles.  Consider 
the  cases  shown  in  Fig.  12,  where  in  a)  the  particles  have  no  transverse 
motion  and  in  b)  they  oscillate  about  the  centroid  of  the  beam.  For  case 
a)  the  spread  in  circulation  frequencies  is  clearly  first  order  in  a/R, 
while  for  b)  the  averaging  over  major  radius  leads  to  a  spread  of  order 
(a/R)2.  For  the  Finn-Manheimer  type  of  equilibrium,23  which  is  the  equi¬ 
librium  we  have  predominantly  used  in  IVORY,  the  beam  behaves  as  in  case  b) 
except  near  the  diamagnetic-paramagnetic  transition,  where  poloidal  rota¬ 
tion  ceases.  The  contrast  between  these  two  cases  can  be  seen  from  the 
IVORY  simulation  results  in  Fig.  13.  In  Fig.  13  (b)  and  (c),  all  para¬ 
meters  are  the  same  except  for  the  energy.  The  diamagnetic-paramagnetic 
transition  energy23  for  the  beam  is  y  =  18.  In  the  simulations,  the  fre¬ 
quency  spread  of  the  particles  shows  up  as  the  spreading  out  of  particle 
bunches  in  the  toroidal  direction,  starting  from  the  initial  positions 
shown  in  Fig.  13  (a).  The  faster  the  bunches  spread  out,  the  higher  the 
circulation  frequency  spread.  For  Fig.  13(b),  y  =  12,  so  that  the  beam  is 
far  from  the  transition,  and  the  poloidal  rotation  leads  to  the  small 
0(a/R)2  frequency  spread  seen  in  the  figure.  (Not  surprisingly,  the  growth 
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be  proven,  of  course.)  Thin  dielectric  linings  or  weak  slow-wave  struc¬ 
tures  in  the  cavity16.17  should  have  similar  effects  but  must  not  signif¬ 
icantly  decrease  the  electromagnetic  cutoff  frequency.  Although  reducing 
growth  for  low  toroidal  modes  tends  to  increase  it  for  high  modes,  beam 
temperature  is  expected  to  control  the  latter.  Implementing  any  of  these 
ideas  does  pose  serious  technical  problems,  however.  Far  off-center  beams 
may  quickly  strike  the  wall,  dielectrics  outgas  and  flash  over,  and  slow- 
wave  structures  are  easily  damaged  by  high  power  beams.  Nonetheless,  addi¬ 
tional  thought  given  to  such  speculative  approaches  may  prove  fruitful. 


Our  theoretical  model  of  negative  mass  instabilities  can  be  extended 
in  at  least  three  attractive  directions.  Properly  treating  periodic, 
strong  focusing  devices,  including  the  stellatron  and  alternating-solenoid 
betatron,  involves  coupling  different  toroidal  modes  through  the  particle 
equations  of  motion.  This  has  already  been  done  for  beam-breakup  instabil¬ 
ities  in  a  racetrack  induction  accelerator,22  and  applying  similar  methods 
to  the  present  model  should  not  prove  difficult.  Secondly,  solving  the 
field  equations  by  boundary  integrals  over  Green's  functions  rather  than  by 
fitting  the  Green's  functions  to  the  boundary,  as  done  here,  may  allow  us 
to  accommodate  more  general  geometries.  Also,  the  massive  cancellations 
observed  numerically  in  several  terms  of  the  dispersion  relation  perhaps 
could  be  performed  analytically  within  an  integral  formulation.  Adding 
temperature  effects  to  the  model  is  the  third,  and  most  difficult,  exten¬ 
sion.  Whether  the  advantages  to  be  gained  by  replacing  the  existing  useful 
but  ad  hoc  methods  for  estimating  stabilization  due  to  temperature  by  a 
more  nearly  self-consistent  description  justifies  the  effort  entailed  is 
uncertain. 
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Figure  10.  Negative  mas 
guide  field 
persion  relal 


We  conclude  this  Section  with  a  brief  comparison  of  dispersion  rela¬ 
tion  curves  with  growth  rates  determined  by  the  computer  simulation  code 
IVORY.  The  simulation  data  points  were,  for  the  most  part,  already  reported 
in  Ref.  6  and  7.  Figure  10  gives  results  for  i.  =  1  for  a  10  kA  beam  in  a  1 
kG  guide  field.  Agreement  is  excellent  except  for  the  y  =  5  data  point. 
The  discrepancy  probably  is  due  to  beam  temperature  in  the  simulation. 
Figure  11  treats  the  same  beam  but  with  y  held  fixed  at  12  and  l  varied 
from  1  to  20.  Again,  agreement  between  theory  and  simulation  is  excel¬ 
lent.  The  few  simulations  performed  with  rigid  disks  instead  of  discrete 
particles  further  vindicate  the  rigid  disk  approximation  made  in  our 
model . 

D.  Summary  and  Discussion 

We  have  developed  an  improved  rigid  disk  model  for  high  current,  low 
temperature  beams  in  conventional  and  modified  betatrons.  The  resulting 
dispersion  relation,  although  quite  complicated,  has  been  solved  analytic¬ 
ally  in  several  useful  limiting  cases  and  numerically  over  a  wide  range  of 
parameters.  A  novel  result  of  this  work  is  the  dominance  of  inductive  over 
electrostatic  fields  at  high  toroidal  mode  numbers,  which  stabilizes  the 
negative  mass  effect  at  high  energies  for  some  toroidal  modes  and  couples 
with  it  for  other  toroidal  modes  to  give  a  hybrid  instability  of  reduced 
growth  rate.  The  limited  utility  of  a  toroidal  magnetic  field  for  slowing 
instability  growth  also  is  noteworthy.  Model  predictions  agree  well  with 
most  available  simulation  data,  but  additional  comparisons  are  desirable. 

Several  possible  techniques  based  on  inductive  coupling  between  the 
beam  and  the  accelerator  cavity  suggest  themselves  for  reducing  or  elimin¬ 
ating  the  negative  mass  instability  at  low  toroidal  mode  numbers  and  moder¬ 
ate  to  high  beam  energies.  We  have  already  seen  that  positioning  the  beam 
near  the  inner  wall  of  the  torus  kills  the  i  =  1  negative  mass  instability 
of  a  10  kA  beam  above  16  MeV.  (That  this  result  for  a  far  off-center  beam 
carries  over  to  beams  and  cavities  with  circular  cross  sections  remains  to 


Figure  8.  Growth  rates  of  the  t  =  20  neoa^ve  mass  instability  versus 
energy  for  10  kA  beams  in  lkG,  lu  and  100  kG  guide  fields. 


satisfied.  For  the  present  parameters  this  is  true  for  Be  less  than  a 
few  kG  at  y  =  10  and  a  few  thousand  kG  at  y  =  50.  Hence,  we  see  a  steady 
drop  in  the  peak  of  the  growth  rate  curve  as  Be  is  increased  above  1  kG 
until  it  obeys  Eq.  (59)  at  100  kG.  Note  that  growth  below  ytr  also  is 
eliminated.  In  contrast,  growth  rates  at  the  highest  energies  shown  in  the 
plot  are  little  reduced. 

Below  yco  instability  behavior  at  A  =  20,  displayed  in  Fig.  8,  is 
qualitatively  similar  to  that  at  A  =  1.  Growth  rates  above  the  inductive 
transition  energy  are  expected,  based  on  Table  1,  to  fall  off  with  magnetic 
field  as  Be-1/2  once  B@  somewhat  exceeds  wo*  This  happens  above  a 
few  kG  for  energies  of  interest  in  Fig.  8.  At  very  high  guide  fields, 
obably  several  thousand  kG,  the  negative  mass  instability  ought  to  be 
eliminated,  leaving  the  inductive  instability  mentioned  earlier  in  this 
Section.  Figure  9  gives  the  growth  rate  at  infinite  field  strength,  for 
which  Eq.  (13)  collapses  to 

Q2  =  e  (69) 

The  instability  sets  in  abruptly  at  yco.  Growth  scales  as  i/R  and,  for 
high  energies,  (u/y3)1/2.  As  a  test  of  our  work,  we  also  solved  the  dis¬ 
persion  relation  in  Appendix  A  with  k  =  n/b  and  obtained  remarkably  similar 
curves. 

Parameters  differing  from  those  in  Table  2  were  run  in  a  few 
instances.  A  10  kA  beam  in  a  1  kG  guide  field  at  l  =  1  and  20  was  consid¬ 
ered.  Changing  the  vertical  field  index  n  from  0.5  to  0.25  or  0.75  modi¬ 
fied  the  overall  magnitudes  of  growth  rates  by  no  more  than  20%,  although 
details  were  somewhat  different  at  low  energies.  Likewise,  increasing  or 
decreasing  the  beam  dimensions  by  a  factor  of  two  caused  no  differences  not 
accountable  in  terms  of  a  rescaling  of  g^,  g^,  and  g^.  Small  changes  in  R 
also  had  negligible  effects.  Moving  the  beam  all  the  way  to  the  inner 
wall,  on  the  other  hand,  eliminated  all  A  =  1  instabilities  above  y  =  34 
but  increased  growth  rates  by  a  factor  of  2.5  at  low  energies.  Growth 
rates  at  i  =  20  were  approximately  doubled. 


40 


.01  )J 


III 


Figure  7. 


Growth  rates  of  the  *  -  1  negative  mass  instability  versus 
energy  for  10  kA  beams  in  1  kG,  10  kG,  and  100  kG  guide 


for  a  cold  beam  go  to  zero  at  y  =  3,  and  there  are  no  unstable  modes  below 
this  energy.  Therefore,  injection  experiments  at  1  MeV  or  less,  and  1  kA 
beam  current  should  not  encounter  the  negative  mass  instability.  This 
assumes  however,  that  the  beam  is  injected  in  one  turn  with  the  parameters 
in  Fig.  14.  The  cutoff  at  y  =  3  is  a  function  of  the  self-fields  of  the 
beam.  At  lower  currents,  this  point  moves  to  smaller  y,  so  that  if  inject¬ 
ion  takes  place  over  many  turns,  the  intermediate  values  of  the  current  may 
be  unstable. 

TABLE  5.  Modified  Betatron  Parameters  in  IVORY  Simulations 


Torus  Major  Radius 

R 

1  m 

Torus  Minor  Radius 

b 

9  cm 

Beam  Radius 

a 

2  cm 

Beam  Current 

I 

1-10  kA 

Beam  Energy 

Y 

3-50 

Toroidal  Field 

Be 

1-5  kG 

External  Field  Index 

1/2 

Upper  Bound  Frequency  Spread 

Ao 

8xl0"5 

The  second  case  we  look  at  is  for  a  10  kA  beam  with  y  *  12  in  a  1  kG 

toroidal  magnetic  field  and  with  Table  5  parameters.  This  is  the  beam  in 

Fig.  13(b)  above.  We  simulated  the  z  =  1  instability,  which  has  a  pre¬ 

dicted  growth  rate  of  5xl0"4  cm-1  for  a  cold  beam.  The  poloidal  rotation 
frequency  is  copr  =  5.7x10" 3  cm-1,  which  exceeds  A0,  so  that  the  beam 
falls  under  case  (b).  The  simulation  results  verify  that  the  small  fre¬ 
quency  spread  seen  in  Fig.  13(b),  has  a  negligible  effect  on  the  growth 

rate.  The  simulation  gives  r  *  5xl0"4  cm-1,  the  cold  beam  growth  rate,  and 
the  be.im  is  strongly  disrupted  in  less  than  1  psec.7 

As  a  third  example,  we  consider  the  beam  in  Fig.  13(c),  with  10  kA  and 
y  =  20.  For  this  case,  the  poloidal  rotation  frequency  is  wpr  =  1.5xl0"3 
cm"1,  which  is  rather  slower  than  for  y  =  12.  However,  for  £  =  1,  the  beam 
still  falls  under  case  (b),  since  a0  <  u>pr.  Thus,  the  simulation  gives 
r  =  2.5xl0"4  cm"1,  compared  to  the  cold-beam  theoretical  growth  rate  of 
2.7xl0-<4  cm"1.  We  note  that  even  if  this  beam  had  the  maximum  possible 


frequency  spread  consistent  with  its  radius,  namely  A0  =  8*10-b  cm-1,  it 
would  still  be  unstable.  In  order  to  make  A0  sufficiently  large  to  kill 
the  instability,  the  beam  radius  must  be  increased. 

We  briefly  mention  an  anomalous  result  encountered  for  this  last  set 
of  parameters.  Since  at  y  =  20  the  beam  is  above  the  transition  energy,  no 
toroidal  magnetic  is  needed  to  confine  it.  Theory  indicates  that  for  Be 
=  0,  the  growth  of  the  instability  should  be  slightly  larger  than  that  for 
Be  =  1  kG.  However,  a  simulation  of  the  Z  =  1  mode  with  Be  =  0  shows  a 
much  smaller  growth  rate,  namely  r  <  3*10-5  cm-1.  The  amount  of  frequency 
spread  seen  is  not  sufficient  to  explain  this  stability.  As  can  be  seen  by 
comparing  Fig.  15  with  Fig.  13(c),  the  frequency  spread  for  Be  =  0  is  in 
fact  significantly  smaller  than  for  B©  =  1  kG.  It  may  be  that  our  theo¬ 
retical  model  is  inaccurate  for  some  reason  when  Be  goes  to  zero.  Since 
the  present  efforts  in  betatron  research  are  concentrating  on  machines  with 
substantial  toroidal  magnetic  fields,  we  have  not  yet  given  much  attention 
to  this  problem.  We  hope  to  be  able  to  return  to  it  in  the  near  future. 
The  examples  we  have  discussed  are  summarized  in  Table  6. 


TABLE  6.  Summary  of  Examples  Discussed  in  Text. 


I  (kA) 

Energy  (  MeV) 

Z 

Stability 

1 

1 

1 

Stable*l  r  =  0  from  Fig.  14 

20 

Stable  ' 

2.5 

1 

Unstable*,  Insufficient  Frequency  Spread 

20 

Stabilized  By  Frequency  Spread 

10 

5.5 

1 

Unstable,  Insufficient  Frequency  Spread 

9.5 

1 

Unstable,  Insufficient  Frequency  Spread 

♦Analytic  prediction  not  verified  with  IVORY. 

The  above  comparisons  between  theory  and  simulations  establish  some 
credibility  for  our  method  of  calculating  the  effective  frequency  spread  on 
a  beam.  We  can  now  compute  upper  bounds  on  the  stable  current,  vmax, 
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Figure  15.  Particle  positions  after  140  ns  for  ■>  -  20  beam  in  zero 
toroidal  magnetic  field.  Compare  with  Fig.  13(c). 


(current  divided  by  17  kA)  once  all  the  other  parameters  are  specified.  We 
do  this  by  assuming  that  the  beam  has  the  maximum  circulation  spread,  A0, 
regardless  of  its  energy.  Results  are  shown  in  Fig.  16.  Note  that  since 
the  growth  varies  as  va,  where  a  <  1/2,  whereas  A0  varies  linearly  with 
beam  radius  a,  a  small  change  in  radius  can  produce  a  large  change  in 

vmax*  Also,  for  small  vmax,  the  current  is  most  limited  at  injection 

(assuming  Yinj  =  3),  because,  as  the  beam  accelerates,  its  stability 
improves.  For  larger  vmaXf  the  beam  may  be  more  unstable  at  energies 
higher  than  the  injection  energy.  This  leads  to  the  crossing  of  the  lines 
in  Fig.  16(b).  This  behavior  is  due  to  the  humped  shape  of  the  growth 
curve,  and  is  related  to  the  multi-valued  nature  of  the  limiting  current  as 
a  function  of  beam  energy  spread  found  by  Sprangle  and  Chernin.5  This 
multi-valuedness  can  occur  even  when  there  is  no  energy  spread.  As  an 
example,  compare  Fig.  16(b)  with  Fig.  14.  Figure  14  shows  that  for  y  <  3, 
a  1  kA  beam  is  stable,  while  Fig.  16(b)  claims  that  I  <  100  A  is  necessary 
for  stability  (assuming  a  2  cm  beam  radius).  The  paradox  arises  from  the 

fact  that  the  growth-rate  cutoff  in  Fig.  14  is  a  function  of  beam  current. 

Thus,  lower  currents  are  unstable  for  y  <  3. 
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bounds  on  stable  currents  in  the  modified  betatron 


IV.  STELLATRON  ACCELERATOR:  EQUILIBRIUM  PROPERTIES 


A.  Particle  Orbits  In  The  Stellatron. 

The  stellatron  accelerator3  adds  strong  focusing  to  the  modified  beta¬ 
tron  device  through  the  use  of  a  helical  quadrupole  field,  as  shown  in 
Fig.  17.  This  results  in  major  improvement  in  the  tolerance  of  the  device 
to  mismatch  between  the  vertical  field  and  the  beam  energy  [see  Refs.  3,  24 
and  Chapter  IVB].  The  work  we  describe  here  arose  out  of  some  IVORY  simu¬ 
lations  of  the  negative  mass  instability  in  this  device.  It  was  observed 
that  for  some  parameters  the  beam  developed  large  envelope  oscillations, 
with  a  few  particles  striking  the  wall  (see  Fig.  18).  This  was  puzzling 
because  previous  results  seemed  to  indicated  that  the  particles  should  be 
confined  for  these  parameters.3  When  we  reduced  the  current  by  orders  of 
magnitude,  we  found  the  same  behavior,  indicating  that  collective  effects 
were  not  involved.  We  then  went  back  to  the  single  particle  equations  in 
the  stellarator  fields,3  and  studied  the  orbits  they  predicted  both 
analytically  and  numerically.  For  the  numerical  work,  we  used  a  single 
particle  orbit  code.  The  details  of  this  investigation  are  given  in 
Appendix  B.  Here  we  summarize  the  main  results. 

The  toroidally  averaged  effect  of  the  quadrupole  fields  is  a  focusing 
one.  At  any  given  toroidal  location,  however,  the  field  is  focusing  in  one 
transverse  direction  and  defocusing  in  the  other.  Because  of  this,  the 
orbit  of  a  particle  can  be  quite  sensitive  to  where  it  starts  out.  Note 
that  the  stability  of  the  orbit  is  independent  of  the  initial  conditions, 
and  so  is  not  in  question  here.  The  effect  of  the  initial  conditions  on 
the  amplitude  of  the  particle  orbit  is  shown  in  Fig.  19.  The  orbits  shown 
are  limiting  cases  in  that  for  one,  the  particle  is  started  in  the  focusing 
plane  and  for  the  other,  the  particle  is  started  in  the  defocusing  plane. 

The  particle  orbits  can  be  computed  analytically  from  the  linearized 
single  particle  equations.  The  general  orbit  has  the  following  character: 
in  a  frame  which  rotates  as  exp(im9/2),  where  m  is  the  mode  number  of  the 
quadrupole  field,  the  orbit  is  the  sum  of  two  right  ellipses  whose 
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Magnetic  field  configuration  in  the  stellatron  device,  showing  the 
superposition  of  a  vertical  mirror  field,  a  toroidal  guide  field, 
and  a  helical  quadrupole  field. 


-7.00  -3.50 
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Figure  19.  Numerically  integrated  orbits  of  particles  which  start  at  a) 
r  =  -2  cm,  2  =0,  9  =o,  b)  z  =  -2  cm,  r=0,  0  =0.  Note 
that  the  dimensions  of  the  boxes  are  different  for  a)  and  b) 


amplitudes  and  phases  are  determined  from  the  initial  conditions.  The 
simplest  initial  conditions  to  analyse  are  those  of  Fig.  19,  where  the  par¬ 
ticles  are  started  off  with  zero  transverse  velocity.  One  can  then  show 
that  for  the  particle  in  the  focusing  (defocusing)  plane,  its  initial  posi¬ 
tion  is  at  an  outer  (inner)  turning  point  of  its  poloidal  motion.  A  plot 
of  the  amplitude  of  the  particle  orbits  versus  the  quadrupole  field 
strength  is  given  in  Fig.  20.  Thus,  the  smallest  amplitude  betatron  osc il¬ 
lations  occur  for  particles  started  in  the  focusing  plane. 

We  note  that  it  is  possible  to  take  advantage  of  the  mismatch  proper¬ 
ties  of  the  stellatron  without  creating  large  envelope  oscillations.  Thus 
Ref.  3  shows  that  for  the  parameters  in  Fig.  18  and  e  =  0.5,  the  beam  can 
easily  withstand  a  50%  energy  mismatch.  Figure  20  shows  that  the  maximum 
excursion  due  to  the  defocusing  of  the  quadrupole  field  is  about  1.6  times 
the  initial  displacement,  for  this  value  of  e. 

One  practical  implication  of  these  results  is  that  a  particle 
injector,  e.g.,  of  the  type  used  in  the  U.C.  Irvine  experiment,  should  be 
located  in  a  focusing  plane  to  ensure  the  smallest  beam  oscillations. 

The  above  calculations  assume  m  >  0,  which  means  that  the  field-lines 
twist  anti-clockwise  as  one  goes  in  the  positive  0  direction  around  the 
torus.  As  pointed  out  in  Ref.  3,  m  <  0  is  an  undesirable  state  from  the 
point  of  view  of  stability.  Even  when  the  orbits  for  m  £  0  are  both 
stable,  however,  m  <  0  gives  larger  amplitude  orbits  than  m  >  0,  as  can  be 
seen  from  Fig.  21.  Thus,  m  <  0  is  also  undesirable  from  the  point  of  view 
of  keeping  the  orbit  amplitudes  small. 

B.  Numerical  Simulations  of  Mismatched  Equilibria  Including  Self- 

Fields. 

In  order  to  investigate  the  self-consistent  behavior  of  a  high  current 
beam  in  a  highly  mismatched  condition,  we  performed  some  IVORY  simula¬ 
tions.  We  chose  a  current  of  1  *A,  Be  =  5  kG,  and  stellarator  windings 
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with  m  =  20  and  e  =  0.7.  The  beam  was  initialized  as  shown  in  Fig.  22. 
Field  components  with  toroidal  mode-numbers  l  -  0,  ±  20  were  kept  in  the 
calculations.  The  nine  groups  of  particles  shown  in  Fig.  22  were  suffi¬ 
cient  to  resolve  the  shortest  wavelength  i.e.,  2nro/20.  The  value  of  the 
vertical  field  was  chosen  such  that  a  beam  with  y  =  3.0  would  be  matched. 
Two  cases  were  looked  at,  one  with  y  =  10,  and  one  with  y  =  1.5.  A  snap¬ 
shot  of  the  beam  with  y  =  10  after  about  seven  revolutions  around  the  torus 
is  shown  in  Fig.  23.  It  can  be  seen  that  even  with  this  huge  mismatch,  the 
beam  is  offset  by  only  about  1  cm,  which  is  consistent  with  the  theoretical 
prediction  obtained  from  Eq.  (7)  of  Ref.  2.  The  latter  predicts  that  if 
the  beam  starts  off  at  r  =  R,  with  y  =  10,  it  should  oscillate  about  the 
position  r  =  R  +  0.7  cm  with  amplitude  of  0.7  cm. 

Figure  20  predicts  that  for  e  =  0.7,  some  particles  should  reach  twice 
their  initial  displacement  from  the  beam  centroid.  In  Fig.  23,  particles 
are  found  at  larger  displacements  than  this,  probably  due  to  the  energy 
mismatch  and  self-field  effects. 

The  results  for  y  =  1.5  after  approximately  two  major  revolutions  are 
shown  in  Fig.  24.  The  beam  offset  due  to  the  mismatch  is  imperceptible. 
However,  some  degree  of  bunching  is  starting  to  occur  in  the  toroidal 
direction  with  wavelength  2uR/m.  This  may  be  due  to  a  stable  longitudinal 
oscillation  of  the  low  energy  ring  excited  by  the  steliarator  windings,  or 
it  could  be  an  instability.  This  simulations  were  not  run  long  enough  to 
distinguish  between  these  two  possibilities. 
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(D)  Negative  Mass  Instability 


As  seen  in  Chap.  V,  the  simulations  we  have  performed  of  the  negative 
mass  instability  in  the  stellatron  show  growth  rates  similar  to  those  in 
the  modified  betatron  at  the  same  parameter  values.  We  have  therefore  com¬ 
puted  the  growth  rates  of  the  negative  mass  instability  for  the  UC I  para¬ 
meters  in  the  absence  of  the  stellatron  windings  in  order  to  give  a  rough 
estimate  of  the  magnitudes  involved.  The  results,  which  assume  zero  fre¬ 
quency  spread,  are  shown  in  Fig.  30.  The  growth  rate  of  the  z  =  1  insta¬ 
bility  is  sufficient  to  make  the  instability  nonlinear  in  1  psec.  We  can 
estimate  an  effective  frequency  spread  on  the  beam  based  on  the  results  in 
Chap.  III.  From  the  single  particle  orbits  in  Fig.  28,  we  compute  a  poloi- 
dal  rotation  frequency  of  u)pr  «  4xl0"3  cm-1.  Assuming  a  2  cm  beam 

radius,  we  have  A0  «  5xl0"4  cm-1.  Thus,  for  £  =  1,  we  have  £A0  <  wpr 

so  that  the  effective  frequency  spread  of  the  particles  is  of  order  (a/R)2 
=  6xl0"5  cm-1.  This  is  insufficient  to  stabilize  the  mode,  since  the  cold 
beam  growth  rate  is  r  •»  2x  10-4  cm-1.  For  £  =  20,  on  the  other  hand,  £A0 
>  u>pp,  so  the  frequency  spread  approaches  the  upper-bound  result  i.e.  A  » 
A 0  =  5x  10-4  cm-1.  The  latter  value  easily  stabilizes  the  z  =  20  mode  for 
which  r  =  l.lxlO"3  cm"1  (from  Fig.  30).  Our  calculations,  therefore,  based 
on  the  assumption  of  a  beam  radius  of  about  2  cm,  predict  that  low 

£-numbers  are  unstable  and  high  £-numbers  are  stable.  Additional  work  is 
necessary  to  confirm  these  predictions,  such  as,  for  example,  an  analytic 
model  which  includes  the  stellatron  fields,  and  IVORY  simulations.  What  is 
clear  from  the  experiment  is  that  if  the  negative  mass  instability  is 
occurring,  it  is  not  causing  any  major  disruption  of  the  beam  since 
acceleration  takes  place  over  ~10Q  ysec.  It  may  be  that  the  strong 

focusing  of  the  stellatron  windings  is  sufficient  to  contain  the  beam  in 
spite  of  the  nonlinear  development  of  the  instability  [cf.  Fig.  25]. 

The  above  calculations  assume  a  monoenergetic  beam.  An  alternative 
explanation  of  beam  stability  is  possible  if  one  assumes  a  significant 
(-  10%)  energy  spread  on  the  beam  and  that  the  beam  radius  is  determined 
solely  by  this  spread.  In  this  case,  there  is  no  net  poloidal  rotation  and 
the  effective  frequency  spread  is  A0,  independent  of  £.  For  a  2  cm 
radius  beam,  this  spread  is  sufficient  to  stabilize  all  £-numbers. 


latron  at  y  =  9  for  energy-matched  particles 


in  Fig.  27.  The  difference  is  greatest  at  low  energies.  Even  at  y  =  9, 
which  is  the  final  energy  of  the  UCI  beam,  the  stellatron  configuration 
tolerates  about  seven  times  more  mismatch  than  the  modified  betatron. 

(B)  Transition  From  Positive  to  Negative  Mass 

The  stellarator  windings  give  rise  to  a  transition  energy  at  which  the 
effective  toroidal  particle  inertia  changes  sign.25  This  is  in  contrast  to 
the  modified  betatron,  where,  in  the  absence  of  self-field  effects,  the 
particle  mass  is  negative  at  all  energies.  The  equation  for  the  transition 
energy  is 


This  equation  is  implicit,  since  y  appears  on  each  side.  To  compute  y$t 
we  rationalized  Eq.  (71)  to  obtain  a  fifth  order  polynomial  in  y.  There  is 
only  one  physically  meaningful  root  in  general,  which  gives  the  transition 
energy.  For  the  UCI  parameters,  yst  =  3.4.  The  beam  is  therefore  in  the 
negative  mass  regime  over  most  of  the  acceleration  cycle. 

(C)  Single-Particle  Orbits 

A  single-particle  code  was  used  to  generate  the  orbits  shown  in  Figs. 
28-29.  The  parameters  used  correspond  to  near  the  end  of  the  acceleration 
cycle.  The  orbits  in  Fig.  28  are  for  a  particle  with  matched  energy.  The 
strong  toroidal  field  keeps  the  particle  confined  to  a  thin  annulus  in  the 
transverse  plane  regardless  of  whether  it  starts  off  in  the  focusing  or 
defocusing  plane.  It  also  gives  rise  to  the  slow  precessional  motion  of 
the  particle  about  the  center  (cf.  Fig.  19,  where  the  particles  precess 
about  the  center  much  faster).  This  slow  precessional  frequency  is  roughly 
one  sixth  of  the  toroidal  frequency.  The  effect  of  10%  and  20%  mismatch  in 
the  energy  is  seen  in  Fig.  29.  These  results  agree  well  with  Eq.  (70)  and 
Fig.  27,  which  predict  radial  shifts  in  the  orbit  of  1.35  cm  and  2.7  cm 
for  the  10%  and  20%  mismatches,  respectively. 


VI.  CALCULATIONS  FOR  THE  UCI  STELLATRON 


In  this  chapter  we  describe  some  calculations  carried  out  specifically 
for  the  stellatron  experiment  at  the  University  of  California  at  Irvine. 
The  most  recent  results  from  this  experiment8  indicate  that  a  beam  of  about 
200  A  is  accelerated  from  the  injection  energy  of  20  kV  up  to  an  energy  of 
4  MeV.  Other  parameters  for  the  machine  are  given  in  Table  7.  For  these 
parameters,  we  have  calculated  the  momentum  compaction,  evaluated  the 
transition  energy  from  positive  to  negative  mass,  looked  at  single  particle 
orbits,  and  generated  some  growth  rate  curves  for  the  negative  mass  insta¬ 
bility. 


Table  7.  UCI  Stellatron  Parameters. 


Torus  Major  Radius 
Torus  Minor  Radius 
Beam  Radius  a 
Beam  Current 
Beam  Energy 
Toroidal  Field 
External  Field  Index 
Stellarator  Fields 
Acceleration  Time 


41  cm 
4  cm 

2  cm  (assumed) 
200  A 

20  kV  -  4  MeV 
0  -  9  kG 
-  1/2 

e  =  0.26,  m  *  12 
100  psec 


(A)  Momentum  Compaction 


We  define  the  momentum  compaction  factor  a  through 


Ar  a  Ay 

R  =  y 


(70) 


where  Ay  is  the  energy  mismatch,  and  Ar  is  the  resulting  radial  shift  of 
the  particle  orbit.  For  the  modified  betatron,  1/a  =  1/2  -  ns,  while  for 
the  stellatron,3  1/a  =  1/2  -  ns  +  p2  (m2  +  mp  -  1/2  +  ns)_1,  where  ns 
is  the  self-field  index,  p  is  the  amplitude  of  the  quadrupole  field  index 
and  p  is  the  ratio  Be/B2.  The  above  definition  of  a  is  for  a  single 
particle.  One  can  define  a  compaction  factor  for  the  beam  as  a  whole  by 
replacing  ns  by  nsa2/b2.  We  compare  the  stellatron  compaction  factors 
to  those  of  the  modified  betatron  configuration,  for  the  UCI  parameters, 
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as  before  at  approximately  the  same  growth-rate.  This  behavior  is  suggest¬ 
ive  of  a  resonant  interaction  between  the  beam  and  a  wave-guide  mode.  The 
analytic  work  performed  for  the  modified  betatron,  however,  indicates  that 
in  that  device,  such  an  interaction  can  never  occur,  because  the  wave-guide 
modes  always  have  higher  frequency  than  the  beam  modes.  The  growth  rate 
predicted  in  the  absence  of  the  helical  windings  (i.e.,  for  the  modified 
betatron)  is  just  r  =  lxlO-3  cm-1,  and  this  value  is  not  sensitive  to  the 
cross  section.  At  this  point,  we  only  speculate  that  one  of  the  transverse 
modes  in  the  stellatron  has  a  frequency  high  enough  to  interact  with  the 
wave-guide  modes  (i.e.,  a  gyrotron-type  resonance). 
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Figure  26.  i  =  20  instability  of  a  10  kA,  *t  =  7  beam  in  a  stellatron 
This  may  be  a  gyrotron-type  instability. 


Nonlinear  development  of  *  =  8  negative  mass  instability 
in  the  presence  of  strong  stellatron  fields  (c  =  1.4). 
Little  radial  particle  loss  is  observed. 


the  absence  of  self-field  effects.25  For  the  parameters  in  the  simula¬ 
tions  just  described,  the  transition  occurs  at  Yst  =  7.9  (see  Chap. 
VI8).  Below  this  energy,  one  would  not  expect  to  see  any  negative  mass 
instability.  (At  high  currents,  other  instabilities  may  exist  below  this 
threshold. ) 

We  now  present  some  results  for  high  current  (~10  kA)  beams  above 
the  negative  mass  transition  energy  and  see  how  the  results  depend  on  the 
strength  of  the  strong  focusing.  We  consider  an  energy  of  5.5  MeV,  with  a 
toroidal  field  of  1  kG,  and  stellatron  windings  with  m  =  8.  The  first  run 
we  perform  is  for  e  =  0.35,  for  which  Yst  =  2.3.  The  simulation  shows 
that  development  of  the  i.  =  8  negative  mass  instability  is  hardly  affected 
by  the  strong  focusing.  The  growth  rate  is  r  «  1.7*10“3  cm"1,  compared  to 
the  modified  betatron  result  of  2*10"3  cm"1.  In  the  nonlinear  stage,  about 
three-quarters  of  the  particles  are  lost  to  the  wall.  When  e  is  increased 
to  1.4,  (yst  =  5*2)  the  results  are  somewhat  different.  The  linear 
1  growth  rate  is  reduced  to  about  7. 5X  10_l+  cm-1,  and  while  the  beam  becomes 
highly  bunched  in  the  toroidal  direction,  little  particle  loss  occurs,  due 
to  the  strength  of  the  focusing  (see  Fig.  25). 

We  report  finally  on  some  stellatron  simulations  which  are  not  clearly 
understood  at  this  time.  These  were  for  a  10  kA,  y  =  7  beam  in  a  5  kG  tor¬ 
oidal  field.  The  stellarator  windings  had  e  =  0.7  and  m  =  20.  In  the 
first  run,  the  torus  was  assumed  to  have  a  square  cross-section,  17.6  x 
17.6  cm.  We  found  that  the  i  -  20  mode  on  the  beam  had  a  large  growth  rate 
r  -  5*10-3  cm-1,  and  the  beam  was  disrupted  in  just  one  major  turn,  as 
shown  in  Fig.  26.  A  subsequent  run  with  a  round  cross  section  with  radius 
b  =  8.8  cm,  and  everything  else  unchanged,  showed  no  growth  in  8  major 
revolutions.  It  was  suspected  that  since  the  1  =  2 0  mode  has  a  frequency 
near  the  electromagnetic  cutoff  the  torus,  wc  =  /2  n/L,  where  L  is  the 
length  of  the  side  in  the  square  cross  section  case,  the  beam  behavior 
might  be  sensitive  to  the  transverse  dimensions  at  this  frequency.  We 
therefore  increased  the  radius  of  the  round  cross  section  case  to  give  it 
the  same  cut-off  as  the  square  cross  section  case,  namely  b  =  2.4  L//2tc  = 
9.5  cm.  With  this  radius,  the  beam  again  exhibited  the  same  strong  growth 
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V.  STELLATRON  ACCELERATOR:  STABILITY 


As  we  saw  in  the  previous  chapter,  the  stellatron  has  an  energy-mis¬ 
match  tolerance  much  superior  to  that  of  the  modified  betatron.  We  now 
attempt  to  discover  whether  this  strong  focusing  has  an  effect  on  beam  sta¬ 
bility.  In  this  work,  we  rely  heavily  on  IVORY  simulations,  lacking  an 
analytic  theory  of  stellatron  stability.  It  must  be  stated  that  the  addi¬ 
tion  of  stellarator  windings  to  the  modified  betatron  adds  two  new  para¬ 
meters,  m  and  n,  to  an  already  large  parameter  space.  A  large  simulation 
code  like  IVORY  is  clearly  not  the  tool  to  use  in  optimizing  the  choice  of 
parameters  in  the  stellatron.  Therefore,  we  confine  ourselves  to  running 
simulations  at  beam  parameters  for  which  we  already  know  and  understand  the 
results  in  the  modified  betatron. 

The  first  simulations  we  describe  are  for  the  Jl  =  20  modes  on  a  mod¬ 
erate-current  (lkA)  beam  at  low  energy  (1-2.5  MeV).  The  other  parameters 
for  these  runs  are  Be  *  5  kG,  R  =  1  m,  b  =  10  cm,  a  =  2  cm,  and  m  =  20 
windings  with  e  =  0.7.  It  is  interesting  to  compare  the  results  with  those 
described  for  the  modified  betatron  in  Chap.  III.  At  1  MeV,  we  find  that 
the  stellatron,  like  the  betatron,  is  stable  through  15  major  revolutions 
(the  length  of  the  simulation).  This  result  is  not  surprising,  since  at 
this  energy  the  cold-beam  dispersion  relation  for  the  modified  betatron 
predicts  stability.  At  2.5  MeV,  we  again  find  stability  in  both  machines, 
while  the  cold  beam  modified  betatron  theory  of  Chap.  II  predicts  instabil¬ 
ity  (Fig.  14).  The  reasons  for  stability  are  different  in  the  two  cases, 
however.  As  argued  in  Chap.  Ill,  the  modified  betatron  run  is  stable 
because  of  the  frequency  spread  on  the  beam.  The  graphical  output  from  the 
stellatron  run,  however,  shows  a  frequency  spread  which  is  almost  an  order 
of  magnitude  smaller.  This  allows  us  to  put  an  upper  bound  of  _<  2.0*10-4 
cm-1  on  the  cold-beam  growth  rate  of  the  instability  in  the  stellatron. 
This  value  is  significantly  less  than  the  growth  rate  predicted  for  the 
modified  betatron  case,  namely  r  =  6xl0”4  cm-1  as  seen  in  Fig.  14.  This 
finding  is  consistent  with  the  fact  that  the  stellatron  has  a  transition 
energy  below  which  the  effective  mass  of  the  particles  is  positive,  even  in 
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VII.  CALCULATIONS  FOR  THE  REVERSING  SOLENOIDAL  LENS  BETATRON 


The  reversing  solenoidal  lens  (RSL)  betatron1*  [Fig.  31],  uses  a  series 
of  alternating-direction  solenoids  to  provide  strong  focusing  of  the  beam. 
An  experimental  machine  is  under  construction  at  the  University  of  New 
Mexico.  Unlike  the  modified  betatron  and  stellatron,  the  average  toroidal 
field  in  this  device  is  zero.  It  is  proposed  to  inject  the  beam  into  the 
toroidal  drift-tube  from  a  shielded  cathode,  i.e.,  so  that  the  beam  has 
zero  canonical  angular  momentum.  This  means  that  when  it  enters  the  mag¬ 
netic  field,  the  beam  rotates  poloidally  at  the  Larmor  frequency  G|_  = 
Ce/2y.  As  the  beam  travels  from  one  lens  to  the  next,  the  direction  of 
rotation  alternates.  Thin  lens  theory26  shows  that  in  order  for  the  orbit 
to  be  stable,  the  "phase-advance  per  cell"  Cl  d/Ve,  where  d  is  the 
length  of  one  solenoid,  must  be  less  than  n.  In  practice,  it  is  advisable 
to  be  well  within  this  constraint.  For  the  experimental  injection  para¬ 
meters,  y  =  2.0,  Be  =  80  Gauss,  and  d  =  30  cm,  we  obtain  Cl  d/Ve  = 
0.4.  Note  that  since  the  poloidal  rotation  of  the  beam  keeps  alternating 
direction,  and  the  phase-advance  is  small,  the  particle  orbits  are  essent¬ 
ially  like  those  of  Fig.  12(a).  The  RSL  therefore  provides  the  maximum 
possible  amount  of  frequency  spread  for  a  given  beam  radius,  which  is  very 
favorable  for  beam  stability. 

The  RSL  betatron  bears  a  resemblance  to  the  bumpy-torus  betatron 
looked  at  by  Chernin  et.  al.27  The  latter  device  consists  of  a  modified 
betatron  with  a  sinusoidal  modulation  of  the  toroidal  field.  As  the  modu¬ 
lation  increases  the  bumpy  torus  becomes  like  a  solenoidal  lens  device  with 
very  thick  lenses.  Stability  calculations  for  the  bumpy  torus  show  that 
the  beam  is  unstable  in  the  vicinity  of  points  where  p2  +  2  -  4  ns  = 
m2q2,  where  q  is  an  integer  and  p  is  the  ratio  of  the  solenoidal  field  to 
the  vertical  field.  Now,  if  the  canonical  poloidal  angular  momentum  of  the 
beam  is  zero,  as  it  is  in  the  solenoidal  lens  device  (since  magnetic  flux 
is  excluded  from  the  cathode),  it  automatical  ly  satisfies  p2  +  2-  4  ns  = 
0,  assuming  that  the  emittance  of  the  beam  is  small.  This  implies  that  the 
beam  lies  at  or  near  the  boundary  of  an  unstable  region  at  all  times,  i.e., 
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Reversing  Solenoidal  Lens  Betatron  configuration.  The  poloidal 
beam  velocity  v  reverses  direction  from  one  solenoid  to  the 


the  beam  is  marginally  stable.  This  is  not  surprising,  since  p2  +  2  -  4 
ns  =  0  is  the  condition  for  Brillouin  flow,  and  therefore  specifies  the 
maximum  current  density  that  can  be  transported  for  a  given  value  of  Be. 
If  Be  drops  below  this  value,  then  the  beam  becomes  unstable.  However, 
all  that  happens  is  that  the  beam  radius  increases  until  the  current  dens¬ 
ity  drops  enough  to  make  the  beam  stable.  If  the  change  in  Be  is  non- 
adiabatic,  then,  in  addition,  envelope  oscillations,  with  frequency  Q|_, 
are  excited.  This  is  what  is  observed  in  IVORY  simulations  of  the  beam  as 
it  travels  around  the  torus.  In  these  simulations,  the  solenoidal  magnetic 
fields  used  are  the  exact,  numerically  computed  values  for  finite-length 
straight  solenoids.  The  solenoids  are  21  cm  in  length,  6  cm  in  radius,  and 
spaced  10.5  cm  apart.  There  are  thus  20  solenoids  around  the  1  m  radius 
torus.  The  value  of  the  solenoidal  magnetic  field  for  which  the  envelope 
oscillations  are  smallest  is  found  to  be  about  100  Gauss  for  a  100  A,  1  MeV 
beam  with  a  2  cm  radius.  As  the  solenoidal  field  is  varied  about  this 
value,  the  envelope  oscillations  increase  in  amplitude.  Runs  at  120  and  50 
Gauss  show  significant  oscillations  (the  beam  envelope  reaches  twice  its 
initial  radius  for  the  latter  case),  but  the  orbits  are  stable  through  one 
major  turn  ( i . e . ,  through  20  lenses),  which  is  as  far  as  we  ran  the  code. 

In  addition  to  simulations  at  the  low  currents  with  which  the  initial 
experiments  will  be  carried  out,  we  have  performed  a  number  of  runs  at  high 
current.  We  find  that  a  solenoidal  field  of  600  Gauss  can  transport  a  10 
kA,  5.5  MeV,  2  cm  radius  beam  with  minimal  (<  10%)  envelope  modulation. 
Note  that  because  of  the  strong-focusing  nature  of  the  solenoidal  lenses, 
the  toroidal  field  contributes  to  the  energy-mismatch  tolerance  of  the 
device.  We  can  compute  the  momentum-compaction  factor  by  noting  that  in 
the  thin-lens,  small  phase-advance  limit,  the  solenoidal  lenses  act  like  a 
continuous  radial  focusing  force  on  the  beam.  The  linearized  radial  force 
equation  is,  in  normalized  units, 

Ar  +  Ar  (1/2  -  ns  +  p2)  V2 z/y2  =  Ay/yR  (72) 
where  Ar  is  the  perturbed  radial  displacement  from  the  matched  orbit,  Ay  is 
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the  perturbed  particle  energy,  and  other  parameters  are  as  defined  in 
Chap.  VI.  Thus,  the  momentum  compaction  factor  is  given  by 

a-1  =  1/2  -  n$  +  p2  (73) 

for  single  particles.  For  the  beam  as  a  whole,  one  replaces  ns  by 
nsa2/b2  to  get  the  momentum-compaction.  Comparing  these  expressions  with 
those  for  the  modified  betatron  in  Sec.  VI,  one  finds  that  for  the  10  kA 
beam,  the  energy  tolerance  in  the  solenoidal  lenses  is  greater  by  about  a 
factor  of  four  than  that  in  a  uniform  toroidal  field. 

Having  obtained  an  equilibrium  for  the  10  kA  beam  with  y  =  12,  we  pro¬ 
ceed  to  perturb  it  slightly,  and  look  at  the  i  =  20  mode  on  the  beam.  We 
find  that  this  mode  grows  rapidly  with  r  =  5.4xl0“3  cm-1.  With  a  growth 
rate  this  large,  the  frequency-spread  damping  £A0  «  1.6xl0"3  cm"1  is 
insufficient  to  give  stability.  The  effect  of  the  instability  on  the  beam 
is  shown  in  Fig.  32.  Classic  negative  mass  (Be  =  0)  behavior  is  exhib¬ 
ited,  with  the  radial  motion  of  the  beam  leading  to  toroidal  bunching.  The 
linear  growth  rate  is,  in  fact,  very  close  to  that  for  a  conventional  (B0 
*  0)  betatron  with  the  same  beam  parameters.  Now,  since  ns  >  1/2  for  the 
beam,  a  self-consistent  equilibrium  does  not  exist  when  B0  =  0  (the  nega¬ 
tive  mass  growth  rate  is  nevertheless  obtainable  from  the  rigid-disk 
model).  Therefore  it  would  appear  from  this  result  that  at  y  =  12,  the 
beam  equilibrium  is  held  together  by  the  solenoidal  lenses,  but  instabili¬ 
ties  behave  as  if  the  toroidal  field  were  zero.  If  this  is  so,  then  the 
advantage  of  the  large  frequency  spread  in  the  RSL  must  be  weighed  against 
the  enhanced  growth  rates  for  Bg  =  0.  These  conclusions  are  speculative 
at  this  point,  and  we  hope  to  perform  additional  analytic  and  numerical 
work  to  clarify  the  linear  behavior  of  the  solenoidal  lens  betatron. 
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APPENDIX  A 


Instability  In  A  Relativistic  Electron  Layer  With  A  Strong 
Azimuthal  Magnetic  Field 
T.  P.  Hughes  and  B.  B.  Godfrey 
Mission  Research  Corporation 
Albuquerque  NM  87106 

Abstract 

A  thin  relativistic  electron  layer  immersed  in  a  strong  azimuthal  mag¬ 
netic  field  between  two  concentric  cylindrical  conductors  is  shown  to  be 
unstable.  The  instability  is  caused  by  the  curvature  of  the  walls,  which 
modifies  the  interaction  between  particles  in  such  a  way  that  the  inductive 
(attractive)  forces  are  greater  than  the  electrostatic  (repulsive)  forces. 
The  growth  rate  of  the  instability  peaks  at  long  axial  wavelengths. 


A  number  of  instabilities  are  known  to  occur  in  relativistic,  rotating 
electron  layers  (Fig.  1),  among  them  the  diocotron ,l  »2  cyclotron  maser,3 
and  negative-mass  4"6  instabilities.  For  these  instabilities  to  occur, 
it  is  essential  for  particles  in  the  layer  to  be  able  to  move  in  the  radial 
direction.  The  instability  discussed  here  differs  from  those  mentioned  in 
that  it  occurs  even  in  the  presence  of  an  infinitely  strong  azimuthal  mag¬ 
netic  field,  which  prevents  any  radial  motion. 

The  idealized  system  we  consider  consists  of  an  infinitely  thin  elec¬ 
tron  layer  of  radius  rb,  rotating  with  velocity  vo  between  two  cylind¬ 
rical  perfectly  conducting  walls,  as  shown  in  Fig.  1.  The  layer  has  infin¬ 
ite  axial  extent,  and  is  immersed  in  an  azimuthal  magnetic  field  which  is 
strong  enough  to  restrict  particles  to  motion  along  the  field  lines.  No 
other  external  fields  are  applied.  The  assumption  that  particles  can  move 
only  along  the  magnetic  field  lines  implies  that  only  those  forces  which 
act  along  the  field  lines  need  to  be  considered  in  the  analysis.  Thus,  the 
layer's  selffields,  and  immobile  neutralizing  ions,  if  present,  have  no 
effect  on  the  dynamics.  If  we  let  the  radius  of  the  layer  become  infinite 
while  keeping  the  wall  separation  ri  -ri  constant,  so  that  the  geometry 
becomes  planar,  then  the  space-charge  oscillations  of  the  electrons  along 
the  field-lines  are  stable.  For  finite  rb,  however,  we  show  in  the  fol¬ 
lowing  that  under  certain  conditions  the  space-charge  oscillations  are 
unstable.  The  curvature  corrections  to  the  field  equations  play  an  essent¬ 
ial  role  in  causing  the  instability. 


a 


While  the  geometry  of  Fig.  1  is  of  interest  from  the  point  of  view  of 
microwave-producing  devices  and  magnetic  insulation,  our  primary  motivation 
in  studying  it  comes  from  present  interest  in  high  current  circular  accel¬ 
erators,  where  a  stronc  magnetic  guide  field  is  used  to  counter  the  space- 
charge  repulsion  of  the  beam  and  to  suppress  negative-mass  instabilities.7 
The  2-D  geometry  of  Fig.  1  provides  a  readily  analyzable  model  of  beam 
dynamics  in  a  strong  circular  guide  field. 

We  first  derive  the  dispersion  relation  for  the  simplest  case,  k=0, 
where  k  is  the  wavenumber  along  the  cylindrical  axis,  and  later  state  the 
result  for  finite  k.  We  use  units  in  which  the  electronic  charge,  e,  and 
mass,  m,  and  the  velocity  of  light,  c,  have  been  scaled  out.  (Velocities 
are  normalized  to  c,  lengths  are  normalized  to  c/u0  =  l  cm,  frequencies 
are  normalized  to  u0,  densities  to  u02m/4ne2,  etc.).  We  use  the  fluid 
equations  for  the  density  n  and  the  azimuthal  velocity  v  to  describe  the 
beam  dynamics. 


3n 


at  +  r  ae  nv 


u 

ao  nv  =  0  , 


(1) 


(L.  +  v  aj 

vat  r  86/ 


TV  =  -  Ec 


(2) 


where  Eq  is  the  azimuthal  electric  field  and  y  =  (1  -v2)1/2. 


-1 
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For  k=0,  transverse-electric  (TE)  and  transverse-magnetic  (TM)  field  com¬ 
ponents  decouple,  and  we  can  write  down  a  field  equation  in  which  only  the 
Ee  field  appears.  Combining  this  equation  with  the  perturbed  azimuthal 
current  density  computed  from  Eqs.  (1)  and  (2),  we  obtain  the  eigenvalue 
equation 


d  L  L  1  ^-1  rE  +  ( 

dr  “2  Zrrj  r  dr  rEe  l 

Lu  -  *  /r  J  \ 


*1  (w  "  AV 


0*) 


E0  =  0  , 


(3) 


where  n0(r),  v0(r)  denote  the  beam  density  and  velocity  profiles,  y0 
=  (1  -  Vo2)-1/2,  and  we  have  assumed  Eq  «  exp(it©  -  io>t).  We  solve 
this  equation  by  assuming  a  density  profile  of  the  form  n0(r)  =  (2v/rt,) 

6(r  -  rb)  where  v  is  Budker's  parameter  (number  of  particles  per  unit 
length  times  the  classical  electron  radius).  Integrating  over  r,  we  obtain 
a  jump  condition  on  dEg/dr.  E@  is  continuous  across  the  layer,  because 
there  is  no  radial  displacement.  Combining  this  with  the  boundary  condi¬ 
tions  on  Eq,  we  obtain  the  following  dispersion  relation. 


(u  -  *vQ/rb)2  +  uv  u2/yoW0(u>)  =  0 

where 


(4) 


we(u.)  « 


0^(wr2)Y^(wri)-Y^(wri)J^(wr2) 


j;(wrb)Y;(Wr1)-Y;(a)rb)j;(wr1)Jj'j;(u)rb)y;(oJr2)-Y;(wrb)j;(oJr2)J 


•i 
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J*,  'll  are  Bessel  functions  of  order  l ,  and  primes  denote  differentia¬ 
tion  with  respect  to  the  argument.  An  approximate  solution  of  t+ns  equa¬ 
tion  is  obtainable  by  expanding  Wq  about  r  =  rt>,  and  using  some  Bessel 

g 

function  identities.  Letting  r^  -  ri  =  a,  r2  -  rfo  =  b,  we  expand 
We  to  second  order  in  toa,  wb,  and  substitute  u  *  lv o/f'b  +  ^  into  the 
resulting  expression.  Keeping  terms  to  order  r2,  we  obtain  the  approximate 
growth  rate. 


In  the  planar  geometry  limit,  r^  -*•  =°,  keeping  a,b,  */rt>,  v/r^ 
fixed,  Eq.  (5)  predicts  stability.  To  get  instability,  the  right  side  must 
be  positive,  which  yields  the  following  criterion. 


The  first  half  of  this  inequality  is  very  similar  to  one  derived  by 
Briggs  and  Neil  as  the  criterion  for  inductive  (as  opposed  to  capacitive) 
coupling  between  the  beam  and  the  outer  wall.6  The  second  half  of  the 

n 

inequality  is  approximate,  because  at  y0  =  ^/(b-a),  the  denominator 
in  Eq.  (5)  vanishes,  due  to  a  resonance  between  the  beam  and  the  lowest  TE 
mode.  The  correct  growth  rate  in  the  vicinity  of  this  resonance,  obtained 
by  keeping  terms  to^rder  r3  in  expanding  Eq.  (4),  is 

r  «  /3  A  (a2b2v/(a+b))1/3/2Yorb2.  (7) 
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As  y02  becomes  larger  than  rb/(b-a),  Eq.  (5)  again  becomes  valid,  and 
the  beam  is  predicted  to  be  stable.  The  criteria  in  Eq.  (6)  imply  that  for 
a  <  b,  the  beam  is  unstable  only  over  a  finite  range  of  Yo»  as  shown  in 
Fig.  2.  The  parameters  in  this  figure  are  rb  =  10  cm,  a  +  b  =  2  cm,  and 
v  =  0.059,  corresponding  to  an  azimuthal  current  of  1  kA  per  cm  (for  v0  ** 
1).  These  values  are  chosen  to  illustrate  the  magnitude  of  the  growth 
rates  one  would  observe  in  a  table-top  size  experiment.  For  a  ^  b,  there 
is  only  a  lower  bound  on  the  unstable  range  of  Yo* 

For  fixed  Yo»  Eq.  (6)  implies  that  the  beam  is  unstable  only  within 
a  finite  region  of  the  space  between  the  walls,  as  shown  in  Fig.  3.  As  one 
gets  close  to  the  inner  or  outer  wall,  the  beam  becomes  stable.  The  growth 
rates  in  Fig.  3  were  obtained  by  numerical  solution  of  Eq.  (4)  and  from  the 
analytic  expressions,  Eqs.  (5)  and  (7).  Equation  (5)  gives  a  good  approxi¬ 
mation  for  the  A  =  10  curve  except  near  the  resonance  point. 

The  first  half  of  the  inequality  in  Eq.  (6)  continues  to  hold  even 
when  the  expansion  in  ua,  u>b  used  to  obtain  it  is  no  longer  valid.  Figure 

3  shows  that  for  l  as  high  as  90,  the  minimum  unstable  value  of  b  is  given 

by  this  inequality.  The  point  of  resonance  (corresponding  to  the  peaks  in 
Fig.  3),  however,  is  given  by  (b-a)/r^  =  1/y02  on^  ^or  1  £  *  <  10.  For 

higher  % ,  the  point  of  resonance  moves  towards  the  outer  wall.  In  Fig.  3, 
we  assume  for  simplicity  that  Yo  remains  constant  as  rb  changes.  In 
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practice,  space-charge  effects  can  cause  y0  to  vary.  Also,  finite  layer¬ 
thickness  effects  are  not  treated  in  our  model.  They  are  expected  to 
become  important  when  the  wavelength  2-n r^/Jl  becomes  comparable  to  the 
layer  thickness. 

For  finite  k,  TE  and  TM  field  components  are  coupled.  An  analysis 
along  the  lines  of  that  previously  described,  but  involving  considerably 
more  algebra,9  leads  to  the  following  dispersion  relation  for  k  *  0 


y  q 
'oH 


(w  -  *Vrb^  +  3^2  roqf  +0^yj 


2  2  2 
k  *  /'b 


2_2 


=  0 


where  q  =  u2  -  k  ,  W©  is  defined  in  Eq.  (4),  and 


(8) 


w  _ Vqr2)Vqrl)~Vqr2)iyqrl) _ 

[JJL^rb>YJt^rl)-YA^rb)JX^rl>j[JAfc’rb)Y7^rz)-YJtfqrb)JJt(qr2) J  * 

As  k  varies  from  0  to  i/r^,  the  growth  rate  drops,  and  the  energy 
threshold  above  which  the  instability  sets  in  increases,  as  shown  in  Fig. 
2.  This  behavior  is  related  to  the  increasing  gap  between  the  beam  mode 
frequency  u>  -  it v0/ r^,  and  the  TE  mode  frequency  u  *  (k2  +  it2 /r^2 ) 1 
as  k  increases.  For  k  =  Jt/rjj,  the  growth  rate  is  zero,  but  a  very  weak 
instability  persists  for  k  >  l/r^. 

Finally,  we  make  the  connection  between  this  work  and  previous  work  on 
the  negative-mass  instability  in  electron  layers.  In  the  notation  of  Ref. 
5,  the  imaginary  part  of  the  frequency  Ei  is  given  by 

C2/c2  =  2Jtv  M/yq  rb2  (b+  +  b_)  ,  (9) 
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where  M  contains  the  dynamics  of  the  layer,  and  b+,  b_  are  the  "wave  admit¬ 
tances".  The  stability  of  the  layer  clearly  depends  on  the  sign  of  both  M 
and  b+  +  b_.  For  M  <  0,  Eq.  (9)  gives  the  usual  negative-mass  instability, 
provided  b+  +  b_  >  0.  The  fact  that  one  can  have  b+  +  b_  <  0,  thereby 
stabilizing  the  layer,  was  first  pointed  out  by  Briggs  and  Neil.6  When  a 
strong  azimuthal  magnetic  field  is  present,  however,  M  >  0,  and  the  condi¬ 
tion  for  stability  is  reversed.  The  layer  is  unstable  when  b+  +  b_  <  0, 
i.e.,  when  Eq.  (6)  is  satisfied.  The  occurrence  of  the  type  of  instability 
we  have  discussed  here  in  a  toroidal  system,  where  the  analysis  is  consid¬ 
erably  more  complex,10  will  be  discussed  in  a  future  publication. 
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FIGURE  CAPTIONS 


Figure  1.  The  thin  layer  of  electrons  of  radius  rt>  and  infinite  axial 

extent  rotates  with  azimuthal  velocity  v0  between  two  concen¬ 
tric  cylindrical  conductors. 

Figure  2.  Growth  rate  r  of  i  =  10  mode  versus  Yo*  Eor  rb  =  9. 95  cm 

(b  >  a)  both  an  upper  and  lower  cutoff  exist.  For  r^  =  10  cm 
(a  =  b)  only  a  lower  cutoff  exists,  r  goes  to  zero  as  k  + 
Jl/rb  (=  1.0  cm-1  for  this  figure). 

Figure  3.  Instability  growth  rate  as  a  function  of  beam  position.  Soli 
and  dashed  lines  give  numerical  and  analytic  results,  respect 
ively.  The  layer  is  stable  near  walls  at  r:  =  9  cm  and  r2  = 
cm.  The  X  gives  the  analytic  growth  rate  at  the  TE  resonance 
(Eq.  (7)). 
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SINGLE-PARTICLE  ORBITS  IN  THE  STELLATRON  ACCELERATOR 
T.  P.  Hughes  and  B.  B.  Godfrey 
Mission  Research  Corporation 
Albuquerque,  NM  87106 


Abstract 


The  stellatron  accelerator  provides  strong  focusing  for  an  intense  electron 
ring  through  the  application  of  a  helical  quadrupole  magnetic  field.  Single¬ 
particle  orbits  in  this  configuration  are  studied.  It  is  found  that  for  certain 
initial  conditions,  the  amplitude  of  the  particle  orbit  increases  as  the  focusing 
strength  is  increased.  Implications  for  particle  injection  into  a  stellatron  are 
discussed. 
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I.  INTRODUCTION 


The  modified  betatron  device  allows  one  to  confine  a  large  current  of 
electrons  using  a  combination  of  a  vertical  mirror  field  Bz  and  a  tor¬ 
oidal  field  .  A  drawback  of  this  device  is  that  the  vertical  field 
must  be  finely  matched  to  the  beam  energy  y  to  keep  the  beam  from  drifting 
off-center  and  striking  the  wall.  The  allowed  mismatch  is  given  by 

^  =  (1/2  -  n$  a2/h2)  b/rQ  (1) 

where  r^  is  the  major  radius,  b  is  the  minor  radius  of  the  torus 
ns  is  the  self-field  index  (ns  =  2v  r02/(a2y3),  where  a  is  the  beam 
radius  and  v  is  Budker's  parameter).  The  stellatron  configuration2*" 
adds  strong  focusing  to  the  modified  betatron  fields  by  applying  a  helical 
quadrupole  magnetic  field  (Fig.  1).  The  effect  of  the  strong  focusing  is 
to  create  a  deeper  and  narrower  "potential  well"  in  the  center  of  the  torus 
than  that  supplied  by  the  weak  focusing  of  the  vertical  field.  The  allowed 
mismatch  in  beam  energy  is  then 

^  =  (1/2  -  ns  a2/b2  +  n2(m2  +  mp  -  1/2  +  ng  a2 /b2 ) -1 )  b/rQ  (?) 

where  p  is  the  ratio  of  the  toroidal  magnetic  field  to  the  vertical  field, 
n  is  the  amplitude  of  the  quadrupole  field  gradient  on  axis,  and  m  is  the 
wavenumber  of  the  quadrupole  field.  By  increasing  the  value  of  u,  the 
allowable  mismatch  can  become  quite  large/*3 
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Figure  1.  Magnetic  field  configuration  in  the  stellatron  device,  showing  the 
superposition  of  a  vertical  mirror  field,  a  toroidal  guide  field, 
and  a  helical  quadrupole  field. 


In  order  to  look  at  the  effect  of  the  stellatron  fields  on  individual 
particles,  rather  than  on  the  beam  centroid,  we  take  the  linearized  equa¬ 
tions  of  motion  for  a  single  particle  and  compute  the  general  solution 
(Sec.  II).  In  Sec.  Ill,  we  obtain  the  amplitude  of  the  orbits  for  some 
specific  conditions,  and  in  Sec.  IV  show  how  the  amplitude  depends  on  the 
strength  of  the  focusing  field.  These  results  are  compared  to  results  from 
a  code  which  numerically  integrates  the  equations  of  motion. 

II.  GENERAL  SOLUTION  OF  LINEAR  PARTICLE  EQUATIONS 


The  stellatron  magnetic  fields  are  those  of  the  modified  betatron, 
namely,  the  vertical  field  Bz  and  toroidal  field  Bt,  together  with  a 
helical  quadrupole  field  of  the  form 

B  =  i  emB^.  (z  cosme  -  r  sinm9) 
r  c  t 

Bz  =  7  EmBt  s’nm9  +  r  cosm6)  (3) 

where  r,  z  denote  deviations  of  the  particle  positions  from  r  =  r0,  z  = 

0,  and  e  =  2n/mp.  From  Ref.  2,  the  linearized  single  particle  equations  of 
motion  in  these  fields  are 


-  R  Z  + 

0 


cosm©  r  +  (iw?  sinm©  z 
o 


0 


z  +  Q  r 
o 


cosm©  z  + 


2 

o 


sinm©  r 


0  (4) 
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rs 


at' 


where  u  =  B  /y,  C  =  B,/y.  Letting  r  +  i  z  «  l  eime/2,  where  m  is  an  even 
o  z  t 

integer,  we  obtain 

£  +  i  (Q  +  mwQ)  t  +  (-^  -  ng  -  5p-  -  +  ^*  =  0  (5) 

„  -iv+t  iv+t 

Next,  suppose  l  -  \  a+e  '  +  b+e  ,  where  a±,  b±  are  complex 

±  “ 

numbers,  then  we  obtain  as  in  Ref.  2, 


~Y~  =  \  -  ns  +  j  p2  +  j  (m  +  p)2  ±  [(y  -  ns  +  j  p2)(m  +  p)2  +  p2]1/2  (6) 

OJ 

0 

The  orbit  is  stable  provided  v±  are  real.  The  coefficients  a+,  b+ 
are  related  by 

2 


pu) 


b.  =  ai 


■r 

1  V+  +  v±  wQ  (p  +  m)  -  </  (y  '  ns  ■  7  mp  '  i 

The  general  solution  can  then  be  written  as 

-i v+t  iv±t 

5  *  1  a+  e  -  +  P+  a*  e 


r~nr 


T 


-j — y-  =  P±  a+ 


(8) 


i  v+t+ 


To  see  what  kind  of  geometric  shape  this  represents,  write  a±  as  a±e 
where  the  <*+  are  real  numbers.  Then  we  obtain 

r  +  iz  =  exp(imw  t/2 )  l  a+  [(p  +  l)cosv+(t-t±)  +  i  (P.,-1 )  si  nv+(t  -t±)j  (9) 

0  +  - 
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Thus,  in  a  frame  rotating  as  exp(iw0t/2),  the  typical  orbit  is  a  sum  of 
two  ellipses  with  major  axes  |a+(P++l)|  and  minor  axes  |a+(P+-l)|.  From 
this  way  of  visualizing  the  orbits,  we  can  at  once  conclude  that  the  ampli¬ 
tude  of  the  orbit  |sj  satisfies 

min{  |  |a+(p++l)|  -  |a_(p_+l)  |  |  ,  |  |a+(p+-l)  |  -  |a_(p.-l)||}  _<  |s| 

1  max{  |o+(p++l)  |  +  |a_(p++l )  | ,  |a+(p+-l)|  +  |a_(p_-l )  | }  (ID) 

III.  INITIAL  CONDITIONS 

The  general  particle  orbit  above  involves  four  unknown  constants  a+ 
and  t+ .  By  specifying  initial  conditions,  we  can  determine  these  con¬ 
stants.  We  look  at  two  simple  cases. 

Case  a).  Suppose  that  at  t  =  0  the  particle  is  a  position  r  =  0 1 * 

z  =  0,  0  =  0,  and  the  transverse  velocities  are  zero.  Then  we  obtain 
t+  =  0,  and 

mu>  in*) 

*  +  •  |50|{v.(p_-l)  +  (P.+DI/D,  O.  =  |50|{v  +  (p+-l)  +  (p++l)}/0 

where 

mw  nui 

o  =  (p++i)[v.(p.-D  +  -£■  (p.+m  -  (p.+D  [V+(p+-d  +  o++Di  (in 
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For  the  particular  Initial  conditions  chosen,  we  can  make  the  following 
ad  hoc  arguments  about  the  amplitude  of  the  orbit.  Since  the  particle  has 
zero  radial  velocity  initially.  It  is  at  one  of  the  radial  turning  points 
of  the  orbit.  (Here  we  are  using  the  term  radial  to  refer  to  |s|). 

Further,  since  the  quadrupole  focusing  force  on  the  particle  is  inward  at 
Its  initial  coordinates,  we  infer  that  the  turning  point  in  question  is  the 
outer  turning  point.  Therefore,  for  these  Initial  conditions,  the  orbit 
amplitude  obeys 

min{ | |a+(p++l) |  -  |o.(p.+l)||,||a+(p+-l)|  -  |a_(p_-l)||}  <  |t|  <  |*0|  (12) 

where  a+  are  given  by  Eq.  (11). 

Case  b).  Now  let  us  consider  a  particle  initially  at  r  =  0,  z  =  |s0| 
6=0,  again  with  zero  transverse  velocities.  The  particle  is  again  at  a 
turning  point,  but,  assuming  the  quadrupole  force  is  stronger  than  the  weak 
focusing  force,  the  initial  net  force  on  the  particle  is  outward.  There¬ 
fore  the  initial  position  is  at  an  inner  turning  point,  and  |S J  obeys 

|50|  1  |s|  £.  max{  |a+(p++l )  j  +  |a_(P_+l)j,  |a  +  (p+-l)|  +  ja_(P_-l)|}  (13) 

where  a±  are  given  by  Eq.  (11)  above  with  the  replacement  ±1  ♦  +1. 


If  we  were  to  start  the  particle  off  at  a  toroidal  position  m6  =  n/2 
instead  of  0  =  0,  then  the  behavior  of  the  two  cases  described  above 
becomes  interchanged.  The  particle  initially  at  r  *  |£0|  z  =  0  obeys 
Eq.  (13),  and  the  one  at  r  =  0,  z  *  |&0|  obeys  Eq.  (12).  This  is  because 
the  quadrupole  fields  have  rotated  90  degrees. 

IV.  RESULTS 


It  can  be  seen  from  Eqs.  (12)  and  (13)  that  the  amplitude  of  the  par¬ 
ticle  motion  is  a  function  of  the  strength  of  the  quadrupole  field,  i.e., 
of  e.  To  see  the  implications  of  this  we  take  some  typical  stellatron  par¬ 
ameters  (r0  =  1  m,  =  5  kG,  y  =  7,  m  =  20),  and  plot  the  turning 
points  obtained  from  these  equations  versus  e  (Fig.  2).  For  case  b),  the 
outer  turning  radius  increases  monotonical ly  with  e,  so  that  the  orbit  can 
intersect  the  wall,  even  though  the  particle  motion  is  "stable"  (v.^  > 

0).  The  numerical  integration  of  Eq.  (4)  for  the  parameters  in  Fig.  2  with 
e  =  1  gives  the  orbits  shown  in  Fig.  3.  The  inner  and  outer  turning  points 
in  this  figure  agree  well  with  the  analytical  results  of  Fig.  2.  Cases  a) 
and  b)  are  limiting  cases.  For  particles  with  the  same  initial  value  of 
|S0|  as  a)  and  b),  but  different  r,z  values,  the  turning  points  lie  with¬ 
in  the  bounds  shown  in  Fig.  2. 

As  one  would  expect,  the  bounds  on  the  particle  orbits  are  linearly 
proportional  to  the  initial  value  of  |S0|.  Apart  from  this,  we  have  not 
tried  to  extract  any  simple  scalings  for  the  amplitudes. 
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Figure  3.  Numerically  integrated  orbits  of  particles  which  start  at  a) 

r  =  -2cm,  z  =  0,  --  =  0,  b)  z  =  -2  cm,  r=0,  c=0.  Note  that 
the  dimensions  of  the  boxes  are  different  for  a)  and  b). 


It  is  clear  from  the  results  presented  that  the  initial  conditions,  as 
well  as  energy  mismatch,  must  be  considered  in  determining  whether  a  par¬ 
ticle  is  confined  or  not.  This  has  practical  implications  for  injection  of 
particles  into  a  stellatron.  In  the  University  of  California  at  Irvine 
experiment,  the  injector  is  located  inside  the  torus,  offset  some  distance 
ri  from  the  center  of  the  minor  cross  section.  In  order  to  create  a  beam 
of  the  smallest  possible  radius,  the  injector  should  be  located  at  a  tor¬ 
oidal  location  corresponding  to  case  a)  above,  i.e.,  where  the  focusing  of 
the  quadrupole  field  in  the  r  direction  is  a  maximum.  In  that  case,  the 
beam  will  have  a  radius  no  larger  than  r* . 

Finally,  it  should  be  noted  that  while  orbit  amplitudes  can  put  an 
upper  bound  on  the  strength  of  the  quadrupole  focusing,  one  can  achieve 
substantial  improvement  in  the  energy  mismatch  characteristics  of  the 
device  without  approaching  this  upper  bound.  Thus,  from  Fig.  2  we  see  that 
for  c  =  0.5,  a  particle  will,  at  worst,  make  excursions  to  1.6  times  its 
initial  amplitude.  As  shown  in  Ref.  2,  however,  this  value  of  e  allows  50% 
energy  mismatch  without  loss  of  confinement. 
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